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ABSTRACT 



The direct problem in transient heat conduction requires the 



determination of conditions at an interior location when conditions 
are known at the boundaries of a solid. Conversely, the inverse 
problem requires the determination of conditions at the boundaries 
of a solid when conditions are known at an interior location. Con- 
sequently special methods are required in the solution of the inverse 
problem. A new method, numerical inversion of the Laplace transform, 
is used to solve this complex problem. Application of this numerical 



technique to the semi-inf init g" cylinder, and sphere is 

made, and the accuracy of solution is discussed. This method of 
solution provides the engineer with a simple, powerful tool that can 
be used in the determination of heat transfer phenomena in a solid. 




Thesis by: Terrill Jay Wendt, entitled "Application of Numerical Inversion 

of the LaPlace Transform to the Inverse Problem in Transient Heat Conduction". 
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I . INTRODUCTION 

Traditionally, heat-transf er problems are concerned with the 

determination of conditions at an interior location when conditions 

are known at the boundaries of a solid. Such problems are known as 

direct problems, and their solutions are obtained through straight- 

* 

forward classical techniques [5] . In many applications (i.e., 
re-entry thermodynamics, rocket motors, quenching), it is desired 
to determine the conditions at the boundaries of a solid when 
conditions are known at an interior location. Problems of this type 
are known as inverse problems and require special methods of solution. 

Many solutions to the inverse problem of heat conduction exist 
in the literature. Stolz [19] uses numerical inversion of a specific 
direct problem to solve for the heat flux at the surface of a solid. 

He notes that there is appreciable damping or lag in time of the heat 
flux from the interior location to the surface. This damping resulted 
in oscillations of the approximate solution with resultant errors. 
Burggraf [4] obtains an "exact" solution to the inverse problem in 
heat conduction with internal heat generation. However Beck [4] in 
commenting on this solution notes that this "exact" method must 
utilize inexact measured temperature data for which accurate time 
derivatives (necessary for this ’’exact*’ solution) are very difficult 
to obtain. In commenting on the inaccuracies involved in using 
measured temperature data, Frank [12] suggests using a least squares 
fit to the data available in order to obtain an essentially smooth 
temperature function for use in the solution of the inverse problem. 
Sparrow, _et. aJ_. [18] use a modification of Laplace transform inversion 

* 

Numbers in brackets refer to references in the bibliography. 
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by applying a convergence factor to the temperature function in the 
s plane. Deverall and Channapragada [11] utilize a convergence factor 
to obtain a solution for the heat flux in an inverse problem. 

It is intended here to present a simple straightforward method 
of obtaining solutions to inverse problems of practical interest. The 
numerical technique involved in inverting the Laplace transform (used 
in solving the inverse problems) was developed by Bellman and Kalaba 
[2 and 3]. This method, based on orthogonal function theory and 
Gaussian quadrature, is developed and tested for accuracy (Tables I 
and II) . This numerical technique is applied to the semi-infinite 
solid, the "long" cylinder and the sphere. The accuracy of solution 
is compared graphically with the known exact solutions. The FORTRAN IV 
programs used in the solutions are listed in Appendix A for use by 
the reader. 
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II. NUMERICAL TECHNIQUE 

INTRODUCTION 

The numerical technique involved was developed by Bellman and 
Kalaba [2] and involves the application of orthogonal function theory 
to the inversion of the Laplace Transform as opposed to using the 
Convolution Theorem which involves integration in the complex plane. 

It is important to note that the method developed in this 
section can be applied to other problems in physics, biology, eco- 
nomics, chemotherapy and radiation [2]. 



THE LAPLACE TRANSFORM 



We will first consider some of the properties of the Laplace 
transform and limitations involved in applying the inverse. 

Let f(t) be any sectionally continuous function defined for 
t >_ 0. Then the Laplace transform of f(t) is defined as 



F(s) = L{f(t)} = 



.00 

e -St f (t)dt 



( 1 ) 



; 0 

where the real part of s is large enough to make the integral (1) 
convergent. 

The Laplacian Operator L is a linear operator and possesses 

many desirable properties, mathematically and computationally. Some 

of these properties are illustrated in the following theorems which 

are included here for completeness: 

'THEOREM I: If L{f (t)} = F (s) and L{f (t)} = F (s) 

11 2 2 

then L{f (t) + f (t)} = F (s) + F (s) . 

12 12 

THEOREM II: If L{f(t)} = F(s) 



then -} = sF(s) - f(0). 

o t 
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THEOREM III: 



If L{f(t)} = F(s) 



then L‘ 



r 3 N f(t)l d N 



3x 



N 



dx 



N 



F(s) 



THEOREM IV: 



If L{f (t) } = F(s) 



then L 



(I 



f(t)dt| = 1 F(s) 






THEOREM V: If L{f(t)} = F(s) 

then L{f(at)} = — F/^-'l. 

a \a / 

THEOREM VI: Initial Value Theorem. 

If L(f(t)} ■ F(s) 
then 

_ sF(s) = f(t). 

g-KO £-XO 

THEOREM VII : Final Value Theorem. 

If L(f(t)} = F(s) 
then 

Limit N Limit . 

_ sF(s) = ^ f(t). 

THEOREM VIII : Superposition Theorem. 

If L{f ( t) } = F(s) and L{g(t)} = G(s) 
then 



(i 



ft 



M f (t)g(t-t)dt| = Lj 



g(t)f (t-t)dt 



« FCs)G(s) 
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INSTABILITY OF THE INVERSE OF THE LAPLACE TRANSFORM 



The instability of the inverse of the Laplace transform can 
be illustrated by considering the function 
f(t) = sin cot. 

The Laplace transform of f(t) is then 



L{sin cot} = 



oo 



( 2 ) 



s 2 +oo 2 

Now it is known that sin cot is bounded in amplitude by (1,-1); but as 
co is increased, sin cot oscillates correspondingly more rapidly. 
However the Laplace transform of sin cot approaches zero as 03 is 
increased. 

Limit 03 



co-k>° 



; 2 +co 2 



= 0 



In other words small changes in s can cause arbitrarily large 
changes in t with respect to the period of oscillation. Consequently 
we cannot sift or filter out rapid oscillations by using values of 
F(s). Such difficulties as mentioned above can be avoided only if 
f(t) is an essentially smooth function, which is an implicit assump- 
tion in the following development. 



NUMERICAL QUADRATURE 



There are several methods that can be used to approximate a 

definite integral in terms of a finite sum. The method employed here 

is based on Gaussian Quadrature and is used because of convergence 

properties and exactness of solution [15]. Consider the integral 

approximation 

rb M 



w(x)f(x)dx = }> A.f(x,). 

J a i=l 1 1 




(3) 

(4) 
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we can assume that the summation involved in equation (3) is exact 
for a prescribed accuracy, e, for any polynomial of degree less than 
or equal to M. The x^ T s and the A^ f s are to be determined so that 
equation (3) is exact for a polynomial f(x) of high degree. Here 
the x^ f s are called the roots (or nodes), the A_/s are called the 
coefficients (or weights) , and w(x) is called the weighting function 
of the function f(x). 

With this background let us examine the numerical approxima- 
tion to an integral which will be used later to invert the Laplace 
transform. Let 

M 

tt(x) = "1 r (x-x ) , (5) 

i=l 

and let P^(x) be a polynomial of degree N. Then we have 

P N (x) = it(x)Q n _ m (x) + P M _ 1 (x). (6) 

Since P^x^ = P M-l^ x i^ 1 = 





rb 


rb 


then 


W (x)?N (x)dx = 

a J 


w(x)P M _ 1 (x)dx + 
a 



w(x)ir(x)Q N _ M (x)dx. (7) 



Now let us choose the x.'s such that 

1 



w(x)TT(x)Q N _ M (x)dx = 0, 



for all polynomials Q^_^(x) . Then from equation (3) we have 



rb 



w 



(x)P N (x)dx = ^ A i P N (x i ). 
i=l 



( 8 ) 



(9) 



If equation (8) is to be valid for all polynomials ^ for 



N M, it must be valid for 
rb 

w(x)7T(x)x K dx = 0, 



K = 1,2,... ,N-M 



Let us now expand it(x) into a polynomial 

TT(x) = C x 1 , 

i=o 



( 10 ) 



( 11 ) 
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i = 0,1,2,. . . ,N-M 



( 12 ) 



so that equation (10) becomes 




which is a system of M equations in the c.'s if N = 2M-1, 

The system of equations in (12) has a unique solution lf_ and 
only if the matrix of coefficients is nonsingular. 

The determinant of this matrix is called a Vandermonde deter- 
minant and can be easily evaluated, since w(x) > 0 in the inverval 
(a ,b) , as 



1 

1 



x . . 
o 



N 



x 

o 



X 



N 

l 



= T7(x-x). (13) 

i>j J 



1 *N ' ' ' *H 

Since the x^ T s are real and distinct points (roots), the determinant 
in (13) does not vanish and the b^'s can be uniquely determined. 

It can be shown that the zeros of tt(x), determined by its 
coefficients (c^) , are real and in the interval (a,b). It can also 
be shown that Q^_^( x ) is a polynomial of degrees M-l, where N = 2M-1. 

Now let tt^(x) and ff\ (x) designate the functions tt^(x) deter- 
mined by their coefficients for M = i and M = j respectively. Then 

f b 

w(x)tt^ (x) 7T , (x)dx = 0 i =/= j (14) 

If i < j, then n\(x) I s a polynomial of degree less than j; and the 
result follows from equation (8), which has been verified. From 
equation (14), TT^x) and 7T^ (x) are said to be orthogonal on the 
interval (a,b) with respect to the weight function w(x). 
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For each value of M, the x^’s in equation (3) are taken as 

the zeros of 7L,(x) as determined by its coefficients. The A u f s in 
M i 

equation (3) can then be determined by the method of undetermined 
coefficients. It can be shown that all of the A^’s are positive. 

Therefore for a given weight function w(x) and interval of 
integration (a,b) , the zeros of the polynomials tt^(x) > as determined 
by their coefficients (c^) , need be computed only once . Likewise 
the corresponding A^ f s need be computed only once . 



GAUSS-LEGENDRE QUADRATURE 



Since we have been dealing with any orthogonal function tt^(x) , 
let us indicate the specific orthogonal functions used in the fol- 
lowing development. The Legendre polynomials are the unique Gauss 
quadrature polynomials [15] in the interval (-1,1) (or any finite 
interval) with weighting function 

w(r) E 1 (15) 

Let us denote the Legendre polynomials, in accordance with 
common practice, by the notation P^(r) where 
P^(r) 7T(r) 2^(N! ) 2 * 

and the interval of expansion is (-1,1). These polynomials are defined 
by the Rodrigue’s formula 



V r) = if h* G 2-1 / • 

From equation (14) it follows that 

f 1 

P^(r)P^(r)dr = 0 M =£ N . 

-l 



(17) 



(18) 



The recursion relations for the Legendre polynomials are then 

P N+I (r) ■ (isi) r V r) - (mi) P N-I (r) • (19) 

and P M (r) - r P„(r) + (^} P„’(r). (20) 
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Consequently the least squares polynomial approximation 
(Gaussian Quadrature) to f(x) over (-1,1) with respect to constant 
weighting function is 



N 

f (r) = Xb K P K (r) ’ 

K=0 

u i_ 2 k+1 

where b = — r — 

1 C z 

-1 



f (r)P (r)dr . 

1C 



Therefore the coefficients (or weights) A, are 
i / \ i 

P N ( r > 

dr , 



A. = 
1 



(r-r ) Pi(r) 
— 1 1 JN 



( 21 ) 

( 22 ) 

(23 



where r_/s are the roots of the Legendre polynomials P^(r). 



SHIFTED LEGENDRE POLYNOMIALS 



We have now examined a specific orthogonal polynomial, the 
Legendre polynomial. However, it is most convenient in the following 
development to consider the shifted Legendre polynomials. The reason 
for employing these polynomials is twofold 0 First, the interval of 
expansion (0,1) of the shifted Legendre polynomials corresponds to 
that involved in obtaining an explicit inverse matrix. Second, the 
coefficients of the shifted Legendre polynomials are all integers , 
which is an important property computationally. 

Since the roots and weights have been tabulated for the 
interval (-1,1) [19], it is a simple matter to obtain the roots and 
weights of the shifted Legendre polynomials in the interval (0,1). 

' Let * 

P N (x) = P N (l-2r) (24) 

* 



where P^ (x) are the shifted Legendre polynomials * 

roots of the shifted Legendre polynomials are 
1+r.j 



x. = — rr 

i 2 



and the weights are 



Consequently the 

(25) 

(26) 
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NUMERICAL INVERSION OF THE LAPLACE TRANSFORM 



We have now laid a proper foundation for obtaining an explicit 
inverse matrix to be used in inverting the Laplace transform. Let us 
consider the approximate inversion of the Laplace transform through 
numerical techniques. 

Let 



F(s) = 



e St u(t)dt. 



(27) 



Now we make a change of the variable of integration to 

-t 



x = e 



therefore 



t = - log x 



Substituting equation (29) into equation (27) we obtain 
x S ^u(-log x)dx = F(s). 



Now let a(x) = u(-log x) 
and then equation (30) becomes 



fi 



x S ^a(x)dx = F(s) 



(28) 

(29) 

(30) 

(31) 

(32) 



s— 1 

If we make g(x) = a(x)x. , then using the quadrature formula 

'1 N * 1 

g(x)dx = A g(x ), (33) 

i=l 



we obtain 



N 



i=l 



A. x . s ^a(x ) = F(s) . 
li i 



(34) 



Letting s = 1,2,...,N, we obtain a linear system of N equations in 



the N unknowns, a(x.), i = 1,2,...,N, 

*L 

> A. x a(x.) = F(k+1) k = 0,1,..., N-l 
i=l 1 1 



It is our intention here to obtain 

N-l 

a ( X -j) = a j|;F ( K+ l) » 
K=0 



i = 1,2,. ..,N, 



(35) 



(36) 



where the a. ? s are constants independent of F(k+ 1) , determined once 
IK 



and for all for each N, and tabulated. 
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( 37 ) 



and 



Let us now make another change of variables such that 

* 

y ± = A ± a(x ± ), 

a K = F(k+1) . 



Equation (35) then becomes 
N 

__ K 

^ X i y i = \ ’ 
1=1 

th 



K = 0,1,. . . ,N-1. 



(38) 



(39) 



Multiplying the K equation by a parameter q , and adding the 

1C 

corresponding terms in all the equations , we get 
_N_ /fr-1 \ N-l 

' 7 " STW 



i=l yK=o 

Now setting 
N-l 

f(x) = ^~q K x " > 
K=0 



K=0 



K 



N 



N-l 



we have 



^, y i £(x i> ■ 

1=1 K =0 



a q , 
k h k ’ 



( 40 ) 

(41) 

(42) 



where f(x) is a polynomial of degree N-l, to be chosen in some 
convenient manner. 

We want to determine the y^’s, and thus we require that 



f(x) = f (x) such that 



f . (x.) = 6 . . 
J i 1 3 



-c 



i + j > 
i = j 



(43) 



an orthogonality condition. If 

f j “ SI Vj * K » (44) 

J K=0 J 

equation (42) reduces to 
N-l 

y j - 51 \q Kj . (45) 

J K=0 J 

where q . , k = 0,1,..., N-l, are determined by the orthogonality 
KJ 

condition in equation (43) . 

Since the interval of integration involved in the preceding 
development is (0,1), we will use the shifted Legendre polynomials 



in approximating the integral in equation (32). Because the x^'s are 
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the zeros of the shifted Legendre polynomial P^ (x) , we can use the 



Lagrange interpolation formula 

V (X) 

f j (X) (x-x.) P N ^( Xj ) 



(46) 



The desired q^ f s in equation (45) are the coefficients in 
this polynomial of degree N-l. Repeating this procedure for each j 
yields the desired inverse matrix (q^). 



INTERLACING OF THE ROOTS 



It should be noted that the roots (x_^) for different N inter- 
lace (Appendix B) . Consequently the values of time 

t ± = -log x ± (47) 
intertwine, and the results for one N may be combined with those of 
another N to obtain both a check on the results and an enlarged set 
of t values for extrapolation purposes. 



CHANGE OF TIME SCALE 



The method of inversion described for the Laplace transform 

provides approximate values for the function u(t) at discrete values 

of t. For general N, we obtain u(t^) at the points t^ = -log x^ where 

* 

x^ are the roots (or zeros) of the shifted Legendre polynomials, P^ (x) . 

It is obvious that the inverse transform is thus determined for 
a very small range of t values. Since increasing N does not signifi- 
cantly increase the range of t, we must resort to another method of 
time scaling. Recall that if 



F(s) = L(u(t)} 



f°° 

e St u(t)dt, 



0 



(48) 



then from Theorem V, it follows that 
L{u(at)} 



(49) 
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( 50 ) 



Therefore using the approximate relation 



N 



i=l 



A i X i S_lu( " a lo § X i } “ 1 F (f) ’ 



we can determine u(t) for values of = -a log x^ e 



ALTERNATE METHODS 



Although the method of inverting the Laplace transform pre- 
sented involves the use of the shifted Legendre polynomials, there 
are several other methods that can be used in obtaining an explicit 
inverse matrix. The first is the Method of Papoulis which utilizes 
the Legendre polynomials and the second is based on a different set 
of orthogonal polynomials, the Laguerre polynomials. 

Papoulis Method 

The Papoulis method of inversion involves the expansion of f(t) 
into a series of Legendre polynomials and uses the classical interval 
of orthogonality (-1,1) • It is clear from this method how a particular 
region of the time scale may be examined by changing the scale factor 
Y [8]. 



Let us make the transformation 



x = e 



y > 0 



where f(t) becomes 



f (t) = f (- — log x) = f (x) , 

Making these substitutions into equation (1) we get 
, fi 



YF(s) = 



x 



<S/Y) ‘ 1 f(x)dx , 



(51) 



(52) 



(53) 



By expressing s in terms of y as s = (2 k+1)y> equation (53) becomes 



YF(2ky+y) = 



fi 



V *A 

x f(x)dx. 



(54) 
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Defining f(x) as an even function (i.e., f (x) = f(-x)), the function 

/N 

f(x) can be expanded into a series of even Legendre polynomials as 
00 

f(x) = ^> C N P 2N (x) * (55) 

N=o 

Using the fact that an even polynomial in e Yt 

of degree 2N, Papoulis shows that the transform of 



the form 



FCs^ = p (s') = ( S ~Y)( S 3y) . . . [s-(2N-l)y] 
^ ' r 2 N'' S '' s(s+2y) ... (s+2Ny) 



(56) 



Replacing s first by y then 3y, 5y, . . . (2N+l)y ... in equation (56), 
gives the following system of equations to evaluate C^: 
yF(y) = C q 
Y F(3y) 

(57) 



2NCi 



2N(2N-2) ... 2C N 



F(2Ny+y) 2 N+1 + (2N+l)(2N+3) + *** + (2N+1) (2N+3) . . . (4N+1) * 
Consequently the C^ f s are determined and the function f (x) in real 



time is 



f(t) = !> C M P 0 ..(e Yt ). 



N=o 



N 2N 



(58) 



Gauss-Laguerre Quadrature 

The Laguerre polynomials are the unique Gauss quadrature 

polynomials in the interval (0,°°) with weighting function 
-ax 



w(x) = e 



(59) 



Denote the Laguerre polynomials by L^(x) where 
L n (x) = tt(x)(-1) N , 



(60) 



and the interval of expansion is (O, 00 )* These polynomials are defined 



by the Rodrigue T s formula 



,N 



T , * ax d 

L N (01X) = 6 d^ 



( X N e -a j 



(61) 
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It follows from equation (14) that the orthogonality condi- 
tion is 

.00 

e ax L^(ax) L^(ax)dx =0 M ^ N (62) 

0 

Consequently the resursion relations are 





l n+i ( x) = (1 + 2N " x) V x) " n 2 l n-i ( x) » 


(63) 


and 


L N+1^ = (N+1)[l'(x) - L n (x)]. 


(64) 



The least squares polynomial approximation to f(x) over (O, 00 ) 



with respect to the weighting function 

, s -ax 
w (x) = e , 



(65) 



is F(x) 

where b 

K 



N 

^ b L (ax) , 



e ax f(x) L (ax)dx . 



( 66 ) 

(67) 



Consequently the coefficients (or weights) become 



-ax 



L N (ax) 



(x-x )L'(ax ) 



dx 



(67) 



0 

Using the properties stated above, a similar procedure to that 
employed with the Legendre polynomials may be performed and an explicit 
inversion formula derived. 



TEST CASES 

A numerical technique has been developed based on shifted 

Legendre polynomials to invert the Laplace transform of a function. 

Let us now examine the accuracy obtained in utilizing this technique. 

Two test functions are considered. The first function is 

u(t) = constant = 1 where the Laplace transform is v(s) = . The 

second test function is u(t) = e t where the Laplace transform is 

v(s) = — 7 — . The results obtained for these test cases are tabulated 
s +1 
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in Tables I and II. Notice that the accuracy decreases as N (number 
of roots) increases. However 8-9 figure accuracy is obtained In 
test cases I and 6-9 figure accuracy is obtained in test case II. 
These results indicate that the numerical inversion technique used 
is very accurate. 
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TABLE I 



rjSJT Cr.SE I 



EXACT VALUE 

0.999999999D 00 
0.999999999D 00 
0. 999999999D 00 
0.999999999D 00 

0 . 999999999D 00 
-.999999999D 00 
0. 999999999D 00 
0.999999999D 00 
0.99^9999990 00 
0. 999999999D 00 
0.999999999D 00 

0.9999999990 00 
0.999999999D 00 
0.999999999D 00 
0.999999999D 00 
0.999999999D 00 
0 . 999999999D 00 
9999999990 00 
0 . 999999999D 00 
0 . 999999999D 00 
0. 999999999D 00 

0.9999999990 00 
0.999999999D jO 
0.999999999D 00 
0 . 999999999D 00 
0.999999999D CO 
0 . 999999999D 00 
0.999999999D 00 
0.9999999990 u 0 
0 . 999999999D 00 
0.999999999D 00 
0.999999999D 00 
0. 999999999D 00 
0.999999999D 00 
0.999999999D 00 



APPROXIMATE VALUE 

0.100000000D 01 
0.999999999D 00 
0.999999999D 00 
0.999999999D 00 

0. 999999999D 00 
0 . 100000000D 01 
0 . 100000000D 01 
0.9999999990 00 
0 . 100000000D 01 
0. 999999999D 00 
0.100000000D 01 

0 . 100000000D 01 
0 . 999999999D 00 
0.1000000000 01 
0 . 100000000D 01 
0.999999999D 00 
0 . IOCoOOOOOD 01 
0.999999999D 00 
0 . 100000000D 01 
0.999999999D 00 
0.1000000000 01 

0.999999966D 00 
0.999999962D 00 
0 . 100000012D 01 
0.100000003D 01 
0.9999999880 00 
0 . 1000000C1D 01 
0.999999891D 00 
0.999999944D 00 
0.100U00004D 01 
0.999999862D 00 
0.1000000090 01 
0.999999977D 00 
0.999999962D 00 
0.100000005D 01 
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TABLE II 

P EOT CA6E II 



EXACT VALUE 


API 


0 . 930568155D 00 


0 . 


0.669990521D 00 


0 . 


0.3300094780 00 


0 . 


0.694318442D-01 


0 . 


0.97455 3956D 00 


0 . 


0.670765 5 9 2D 00 


0 . 


0.702922575D 00 


0 . 


0.499999999D 00 


0 . 


c . 297077424D 00 


0 . 


0.1292344070 00 


0 . 


0 . 25 4 4 o04 380-01 


0 . 


0.96o953264D 00 


0 . 


0.9325316830 00 


0 . 


0.6397047840 00 


0 . 


0.716697697D 00 


0 . 


0.6744371690 00 


0 . 


0.425 ^626300 00 


0 . 


0.2633023020 00 


0 . 


0 . 160295 215D 00 


0 . 


0.674683166D-01 


0 . 


0.1304o7357D-01 


0 . 


0.9931419040 00 


0 . 


0. '.642174410 00 


0 . 


0.9136006570 Ow 


0 . 


0.843646452D 00 


0 . 


0.757624318D 00 


0 . 


0.6595561840 00 


0 . 


0.5540274740 -.0 


0 . 


0 .4459725 A5D 00 


0 . 


0.3404438150 00 


0 . 


0.242375681D 00 


0 . 


0.1563535470 00 


0 . 


0.663993424D-01 


0 . 


0.3576255810-01 


0 . 


0.6656095650-02 


0 . 



‘HOXIMA i?iS VALUE 

930568155D 00 
669990521D 00 
3300094780 00 
6943184420-01 

974553956D 00 
870765592D 00 
702922575D 00 
5000000000 00 
2970774240 00 
1292344070 00 
2544604380-61 

98695 32640 00 
9325316830 00 
8397047840 00 
7166976970 00 

574437I69D 00 

425562830D 00 
2833023-./2D 00 
1602952150 00 
6746831630-01 
130467 358D-01 

9931418330 00 
964217450D 00 
9136007610 00 
643646317D 00 
7 5 76 2 4 35 0D 00 
659556165D 00 
5540273690 „0 
4459725870 00 
3404-37930 00 
2423756490 00 
1563535510 00 
863993391D-01 
3578256070-01 
6858151400-02 



III. THE INVERSE PROBLEM IN TRANSIENT HEAT CONDUCTION 



INTRODUCTION 

The problems presented in this section are representative of 
typical one, two, and three dimensional geometries. They were selected 
because of the existence of exact solutions, necessary for comparison 
of results. However, more complicated geometries can be accommodated 
by this method of solution, but the equations and analysis become 
quite complex. The FORTRAN IV programs used in the solution of these 
problems are listed in Appendix A. 



DIMENSIONLESS VARIABLES 

The solution to a problem in heat conduction can always be 
expressed in terms of dimensionless variables. This is usually done 
when making numerical calculations and is adopted here because of 
simplification in the differential equations involved and generality 
in application of the solutions obtained. 

The dimensionless variables used in the following problems are 
Kt Kt 



where 



T = 



K = 



P = 



u = 



£2 R2 

k 

Cp (density) 5 

x _ r 

i " R 

uK uK 
$R 2 



tt N R Y N 
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THE SEMI-INFINITE SOLID 



The method of numerical inversion of the Laplace transform 
is now applied to the solution of the semi-infinite solid. Since 
there exists no physical characteristic length, dimensional variables 
are used in the solution to this problem. 

The partial differential equation governing the temperature 
function is 



9u 9 z u 

3t K Ite 7 ’ 




(69) 


with the boundary conditions 






u(0,t) = Bt , 




(70) 


and 






u(x,0) = 0 




(71) 


The solution of this 


direct problem is 




u(x,t) - 43ti 2 erfc - 


X 


(72) 



JTk7 

Using this solution as the known temperature function at x^ (an 
interior location), we can solve the inverse problem to obtain u(0,t) 
at x = 0 (the surface of the solid). Here the exact solution is known 
as 



u(0,t) = 3t . (73) 

Transforming the differential equation into the s plane we obtain the 



ordinary differential equation 

d 2 v s n 
d?-K V = ° 



(74) 



where v is the Laplace transform of u(x,t) and is a function of x and 



s. The solution to equation (74) is of the form 
v(x,s) = A e x 

Solving this equation at x = x^ we find that 

, N Js/K x, 

A “ v(x 1 ,s)e' r 1 



(75) 



(76) 



3Q 



( 77 ) 



and therefore the solution of equation (74) becomes 

, x f x “ vJs/K (x-x ) 
v(x,s) = v(x^,s)e N 1 

Since the temperature function at an interior location x^ is given 

in equation (72), we can solve equation (77) by taking the Laplace 

transform of equation (72) which then becomes 

vCx^s) = “2 e X 1 (78) 

Substituting equation (78) into equation (77) we have the solution 

in the s plane as 

r( v .) = (§* sP* (x - x l ) (79) 



V V 

At the surface of the solid x = 0, and the solution is 

x. 



V 



<0,s) -/§, e-J^ x l)e^ 



(80) 



Taking the Laplace inverse of equation (80) we obtain the solution 
in the time plane at x = 0 (the surface of the solid) . 

An illustration of the semi-infinite problem is presented 
in Figure 1. In obtaining the results plotted in Figure 2 it was 
found advantageous to use a convergence factor. The convergence 
factor e t was added to the function u(x,t) to improve the accuracy 
of solution. The function u(x,t) was multiplied by e ^ transformed 
into the s plane and subsequently inverted back into the time plane. 
The final answer was then multiplied by e t . Note that the Laplace 
transform of this function is 

L{e t u(x,t)} = F(s+1). 



THE LONG CYLINDER 

The basic assumption involved in the solution of the ’’long" 
cylinder is that the cylinder is sufficiently long so that end effects 
can be neglected. In other words conduction is considered to be only 
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in the radial direction (Figure 1) . The governing partial differential 
equation for the temperature function u(r,t) is 



_9u 

3t 






subject to the boundary conditions 
u(R,t) = 6t , 

and 



(81) 

(82) 

(83) 



u(r,0) = 0 

In terms of dimensionless variables, the solution of the direct 
problem can be represented as the sum of an infinite series as 

J o(paN) 



u(p,T) = T - + 2 *> e 

N=1 

where the are defined by 

J 0 (a N > = 0 • 



-a N 2 x 



a N 3 J,(a N ) 



(84) 



(85) 



Using this solution as the known temperature function at 
(an interior location), we can solve the inverse problem to obtain 
u(p,T) at p = 1 (the surface). Here the exact solution is known as 
u(l,T) = T . (86) 

Proceeding as before we find that the solution of the inverse 
problem in the s plane is 



v(l,s) = v(p 



where 



i’ s) ■ 



v(p, ,s) 



= ±_khl + 



f 1 \ Jq(P 

i^N 2+s / a N d J 1 



(Pj a N) 
(^n) 



(87) 

( 88 ) 



1* s 2 4s 

By numerically inverting equation (87) , we obtain the solution of 
the inverse problem at the surface (p = 1) where the exact solution 
is given by equation (86) . 

The results obtained for the "long" cylinder are plotted in 
Figures 5 and 6. It should be noted that as the distance from the 
surface is increased, the accuracy of solution decreases (Figure 5). 
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This is reasonable since the temperature will be dampened more as the 
distance is increased. Let us consider a characteristic parameter 



then physical significance can be given to V}« Now we can examine the 
accuracy of solution as r) is varied. In Figure 5, surface temperature 
is plotted versus ri for 6 = constant. It can be seen that as ri is 
increased, the accuracy of solution decreases- Since the only variable 
in r) is now r, it is obvious that as r increases, the accuracy of 
solution increases. By interpreting small n as "close" to the surface 
and large ri as "far away" from the surface, we can see the effects 
of damping on the solution. 



it can be seen from Figure 4 that the solution is improved. However 
the results obtained are not quite as good as the other problems 
considered. The inaccuracy introduced into this problem is due to 
the evaluation of an infinite series of Bessel functions. 

THE SPHERE 

The basic assumption involved in the solution of the sphere 
is axial symmetry. In other words the temperature function is not 
dependent on <j) (Figure 1) . The governing partial differential equa- 
tion for the temperature function u(r,t) is 



1-r 



n 




If a thermal penetration depth is taken as 



6 = 2jKt , 



A convergence factor was also used in this problem to improve 
the accuracy. As before, e C was used as the convergence factor, and 




( 89 ) 
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subject to the boundary conditions 

u(R,t) = 8t , (90) 

and 

u(r,0) = 0 . (91) 

In terms of dimensionless variables, the solution of the 
direct problem can be represented as the sum of an infinite series 

as M 

- , N 1-p 2 2 ^ (-1) N -N 2 tt 2 t x /MX 

u (p,T) = T - -f- - -IT vr 3 e sin(NTTp) . (92) 

P N=1 

Using this solution as the known temperature function at 
(an interior location) , we can solve the inverse problem to 
obtain u(p,T) at p = 1 (the surface). Here the exact solution is 
known as 



u (1 ,t) = T . (93) 

Proceeding as before we find that the solution of the inverse 
problem in the s plane is 

pi sinh (Js) 

v(l,s> - v(p r s) sinh Pi> . (94) 



where 



, s 1 l-Pl 2 2 T7" (-1) H ( 1 \ . w „ tN 

v(Pi,s) - - -jj- - ^ 2. -p— (p^) slI '( 1, "Pl>-<95) 



By numerically inverting equation (94) , we obtain the solution of 
the inverse problem at the surface (p = 1) where the exact solution 
is given by equation (93) . 

The results obtained for the sphere are plotted in Figure 7. 
Here as before a convergence factor e ^ was used to improve the 
solution. The results shown in Figure 7 are in very good agreement 
with the known exact solution. 
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INVERSE PROBLEMS IN TRANSIENT 
HEAT CONDUCTION 



SEMI-INFINITE SOLID 

ASSUMPTION 

CONDUCTION OCCURS 
ONLY IN X-DIRECTION 





TEMPERATURE FUNCTION DOES 
NOT DEPEND ON <p. THERE- 
FORE CONDUCTION OCCURS 
IN THE RADIAL DIRECTION. 




NEGLECTED AND THUS 

CONDUCTION OCCURS ONLY 
IN THE RADIAL DIRECTION. 



FIGURE I 
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( 0, T ) (DIMENSIONLESS SURFACE TEMPERATURE.) 



SURFACE TEMPERATURE vs. TIME 

FOR 

SEMI-INFINITE SOLID. 




FIGURE 2. 
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(0,T ) (DIMENSIONLESS SURFACE TEMPERATURE.) 



o 
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SURFACE TEMPERATURE vs. TIME 

FOR 

"LONG" CYLINDER 



A 



□ 



13 




FIGURE 3. 
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u ( 0,T ) (DIMENSIONLESS SURFACE TEMPERATURE.) 



SURFACE TEMPERATURE vs. TIME 
FOR 

"LONG" CYLINDER 
( CONVERGENCE FACTOR ) 




FIGURE 4. 
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SURFACE TEMPERATURE vs. n FOR "LONG" CYLINDER 

(CONVERGENCE FACTOR AND K= I) 
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FIGURE 



u (0,r ) (DIMENSIONLESS SURFACE TEMPERATURE.) 



SURFACE TEMPERATURE vs. TIME 
FOR 

SPHERE 




T ( DIMENISONLESS TIME ) 
FIGURE 6. 
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DISCUSSION OF RESULTS 



In the determination of surface temperature by the method of 
numerical inversion of the Laplace transform, it was found necessary 
to utilize a convergence factor to improve accuracy* The function 
u(x,t) was multiplied by e t , transformed into the s plane, and 
numerically inverted into the time plane. Then the function obtained 
was multiplied by to obtain the desired surface temperature function 
u(0,t)* This method caused the approximation to converge to the known 
exact solution. In the semi-infinite solid and the sphere, the 
results obtained from the use of the convergence factor were very 
good. In the "long" cylinder the results obtained were not quite as 
good. The inaccuracy introduced in this problem is due to the evalua- 
tion of an infinite series of Bessel functions. 

The damping effects were considered for the case of the "long" 
cylinder. It was found that as the distance from the surface was 
increased, the accuracy of the solution decreased. This was a 
reasonable result, since it was logical to expect the solution to 
lag as the distance from the surface was increased. 

The use of time scaling (mentioned in section II) was only 
incorporated in the semi- inf inite solid FORTRAN IV program, but it 
can be included in the other problems. Although the initial tempera- 
ture distribution was taken as zero, any reasonable initial tempera- 
ture distribution can be handled with this method. Such a change 
would simply result in more terms to be inverted. Also it should be 
noted that although the solution was obtained for the surface tempera- 
ture, heat flux at the surface can also be determined. 
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CONCLUSIONS 



Three geometries were considered in the solution to the 
inverse problem in transient heat conduction. These geometries, 
the semi- inf inite solid, "long" cylinder, and sphere, are typical 
of one, two, and three dimensional problems respectively (it should 
be noted here that an analogous approach can be made in the solution 
to the finite plane bounded by two parallel planes) . Solutions were 
obtained which agreed very well to the known exact solution. 

Several parameters were investigated with respect to the 
effects on the accuracy of solution. It was found in all cases that 
utilizing a convergence factor caused the solution to better approxi- 
mate the known exact solution. The damping effect was considered in 
the "long" cylinder problem. As expected it was found that as 
increased (distance from the surface decreased) the accuracy of 
solution increased. 

Therefore the solutions obtained indicate that the method of 
numerically inverting the Laplace transform provides an accurate 
approximation to the known exact solutions. 
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IV. SUMMARY 



The numerical technique involved in numerically inverting the 
Laplace transform was methodically developed. Elementary theorems 
of the Laplace transform were introduced, and the instability of the 
inverse of the Laplace transform was briefly discussed. Gaussian 
quadrature methods were subsequently introduced and developed. The 
shifted Legendre polynomials were then examined for several reasons. 

It was found that the coefficients of the shifted polynomials were 
integers, and the region of orthogonality (0,1) was found to be most 
convenient in the following development. Then, an explicit inverse 
matrix was found which would be used later to solve the inverse prob- 
lem in transient heat conduction. Two useful numerical properties 
were discussed. Interlacing of the roots (and therefore times) was 
found useful in interpolation of the results, while change of the 
time scale was found useful in extending the range of times. Two 
alternate methods of solution (Papoulis and Gauss-Laguerre quadrature) 
that could be used to obtain an explicit inverse matrix were subse- 
quently introduced. Finally, the explicit inverse matrix was tested 
for accuracy. The two test cases employed were for the functions 

_ t 

u(t) = 1, and u(t) = e . The accuracy was found to be very good for 
both cases. 

Three geometries were considered in the solution of the 
inverse' problem in transient heat conduction. These problems, the 
semi-infinite solid, "long" cylinder and sphere, are typical of one, 
two, and three dimensional problems * They were chosen because of the 
existence of exact solutions. It should be noted that the results 
obtained were for discrete points and not for a continuous spectra of 
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information. Consequently this method is applicable to problems 
where discrete information is known (i.e., data from experiment). 

Several parameters were examined with respect to influence 
on the accuracy of solution. A convergence factor e t was found 
to improve the accuracy of solution in all cases considered. It 
was found that as the distance from the surface was increased, 
the accuracy of solution decreased. This was considered reasonable 
since damping effects should be more noticeable at greater distances 
from the surface. 

Consequently the solution to the inverse problem in transient 
heat conduction by numerical inversion of the Laplace transform was 
found to yield accurate results for semi- inf inite solid, "long" 
cylinder , and sphere . 
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APPENDIX A 



FORTRAN IV PROGRAMS 

1. The Semi- inf inite Solid 
The following program is listed: 

SEMI-INFINITE PLANE SOLUTION USING DIMENSIONLESS 
VARIABLES IN ONE-DIMENSIONAL SYSTEM. 

The following subroutines are needed: 



INPUT 


Subroutine 


EXPINV 


Subroutine 



47 



FORTRAN IV PROGRAM SYMBOLS 



MAIN PROGRAMS 

A Explicit inverse matrix 

R Roots of the shifted Legendre polynomials (x„) 

* 

W Weights of the shifted Legendre polynomials (A_^) 

T1 Values of time (t e = -log x.) 

1 & i 

Kt Kt 

TAU Non-dimensional time (t = -jj - 

RO Distance (p) 

R01 Interior location (p^) 

APPR Approximate value of v(x,s) 

XK Reciprocal of K 

EXACT Known exact solution 

SCALE Scale factor (a) 

C Value of e t u(p,T) 

ALFA Roots of the equation J q ( ct^) = 0 

ARG Argument of the function d 0 (^p]_) 

ARQ Argument of the function p^) 

ARP Argument of the function I ( ) 

BESSO Resultant value of 

BESS1 Resultant value of 

BISS Resultant value of I ( Ss ) 

o ^ 

BISS1 Resultant value of I (Js p-) 

o ^ 1 

ADA Value of r) 

N Number of roots 
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EXPINV SUBROUTINE 



C Legendre coefficients 

B Coefficients of the polynomials 

P Legendre derivatives 



, N (x)/(x-x i > 
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C SOLUTION OF THE INVERSE PROBLEM IN TRANSIENT HEAT CONDUCTION 
C 

C 60UN0ARY CONDITIONS ARE: 

- C 1. INITIAL TEMPERATURE OISTRIRUTION IS ZERO. 

C 2. SURFACE TEMPERATURE DISTRIBUTION IS: U ( 1 , T AU ) =T AU 

C 

C I. SEMI-INFINITE PLANE SOLUTION USING DIMENSIONLESS VARIABLES 

— C IN-ONE-OI MENS IONAL SYSTEM. 

C 

DIMENSION Tl(15).A(15,l5),R(15),W(15)*TAU(15),C(15),RO(ll),APP 
*R< 15) , EXACT! 15 > ,F( 20) 

R E AL * 8 A , R , W , C-»-T 1-t-TAlFf-ARPfl-rEXA GT , F, XK, RO. AOA <^0 ) ♦^r-SGAL-F-rT^rtTM 1 — 

*R(20> » R01 (201 

DO 100 J J = 2 * 15 
CALL INPUT(N,R,W,T1) 

CALL EXPINV(fwTWW-rA-) 

WRI TE( 6, 106) 

DO 20 1=1, N 

20 WRITE(6»102)R( I ) ♦ W ( I ) 

SCALE= 2.00*00 

XK= 1 . 0D + C2 
PI=3. 1 A159265 
DO A I 1 = 1, N 

A T AU ( I I4-=SC AL E* I-1L444/-XK 

IP=1 

R01( IP)=C. 00*00 
I R = 1 

R0( IR)=0,OD*00 

2A CONTINUE 

DO AO I S= 1 , N 
Q=IS+1 

S=Q/ SCALE 

EPPR ( I SI=DEXP(-DSQRT(S*XK)*R01< IP) )/S**2 

APPR ( IS)=EPPR(IS)*DEXP<-DSQRT<S*XK)*(RO( IR)-R01< IP) ) )/<XK*SCAL 

*E) 

EXACT< IS)=TAU< IS) — 

ADA( I S ) = R01 1 IP)/<2.0D+00*0SQRT(T1(IS)/XK)) 

AO CONTINUE 

DO 61 I = 1 , N 

C ( I ) =0.CD*00 

DO 61 J= 1 , N 

61 C ( I > = C< I )*A( I, J)*APPR< J) 

DO 62 J=1,N 

62 F ( J ) =C ( J ) *DE XR-m <^)44 

WRI TE! 6, 83 )N 

WR I TE ( 6 , 80 ) RO 1 ( IP) 

WRITE< 6,60 ) 

WR I TE <6,-70 ) I-T-H-I-)-t-F-<-I-)-, E-X A C T ( I ) , ADA( I ) , I ^W-N-) 

IP=IP+i 

R01 ( I P ) = R01 ( IP-D+O. 1D*00 
IF (ROl(IP).LT. 1.000) GO TO 2A 

100 CGNTINUE 

CALL EXIT 

60 F0RMAT(//,11X, 'TIME', IAX, 'APPROXIMATE VALUE '♦ 1 5X ,' EXACT VALUE' 

*,20X, 'ADA' ,//> 

7 0 F 0 R M A T 1 4 < 2 X ,-024 . 1-7-) , /-) 

80 FORMAT ( //, 3X, 'THE VALUE OF R01 IS ',F6.4,//) 

83 F0RMAT(//,3X, 'THE VALUE OF N IS ',15,//) 

102 FORMAT!//, 2<2X,02A. 17) ,/) 

106 FORMAT ( //. 10X. 'ROOTS 13X, 'WEIGHTS'-,//) 

108 FORMAT ( //» 5X, ' T =-LOG R'* 20 X,/, 02 A. 17,/) 

ENO 
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2. The "Long" Cylinder 



The following program is listed: 



TWO-DIMENSIONAL "LONG" 


CYLINDER SOLUTION 


The following subroutines 


are needed: 


INPUT 


subroutine 


EXPINV 


subroutine 


BESI 


subroutine 


BESJ 


subroutine 
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2 



20 



4 



24 



6 



90 



5C 

62 



100 

60 

70 
80 
83 
1C 2 
106 
1C8 
110 



SOLUTION OF THE INVFRSE PROBLEM IN TRANSIENT HEAT CONDUCTION 

BOUNDARY CONDITIONS ARE: 

1. INITIAL TEMP DISTRIBUTION IS ZERO. 

2. SURFACE TFMP DISTRIBUTION IS : U ( 1 , T AU ) =T AU . 
III. TWO-DIMENSIONAL "LONG*'C YL I NDER SOLUTION. 



OIMFNSION T 1 ( IS) , A ( 18.15). R( 15I.WI 15) , T AU ( 1 5 ) ,C ( 18 1 .RD( IS). APP 
*R( 15) .EXACT! 15 ), ARG! 20) , SUM! 15) ♦ ALFA ( 15 ) , ARP( 20 ) , A RQ ( 2 0 » , RO 1 ( ?0 ) 

PFAL*R A.R.W.C.T1 ,TAU , APPR , EX ACT , ALFA . SUM , BE SSO. BESS 1 , B I S S , B I S 
*S1,ARG,ARQ,ARP,S,XK,F(20) 

DO 2 1 = 1, 10 - - 

R F AO ( 5 « 1 10)ALFA( I I 
DO 100 J J = 2 » 15 
CALL INPUT (N,R»W»T1) 

CALL E XP I NV !N»R,W»A) - - 

WRITE! 6»1 06) 

DO 20 I = 1 » N 

WRI TE( 6. 102 )R( I), HU) 

XK=1. 60 + 01 - 

DO 4 I 1 = 1, N 

TAUC I I) = T1< I I )/XK 
PI=3. 14159265 

IR = 1 - • 

R 0 ( I R ) = 1 .00*00 
I P = 1 

R 0 1 ( I P ) =0 .00+00 

CONTINUE 

no 90 i s= i , n 

SUM! I S ) =0.00+00 



S=I S+l 

DO 6 K = 1 » 10 • - 

ARG ( K ) =R01 C IP)*ALFA!K) 

CALL BESJ! ARG, 0,B ESSO, l.OD- 12.1 ER) 

CALL BESJ (ALFA, l.BESSl. 1,00-12, I ER) 

SUM! I S)=SUM( IS ) + < 1.00+00/! ALFA(K)**2+S) )*BESSO/< ALFA! K )**3*BES 

*S1 ) 

ARP! I S ) =DSORT ( S ) 

ARCH I S >=R01 ( I P )*DSQRT (S ) 

- CALL BESI! ARO.C.BISSl, IER) - 

CALL BESI (ARP, 0, BISS, IER) 

APPR ( I S) =( ! 1.0D+00/S**2 )-( ! 1.0D+00-R01 ( I P )**2 ) / ( 4.00+00*S ) > +2. 
*00+00* SUM! IS) )*BISS/(BI$SI*XK) 

EXACT! IS) = TAU( IS) ' 

CONTINUE 
00 50 1=1, N 
C! I )=0. 00+00 

- DO 50 J=1,N - - - 

C c I )=C( I ) + A ( I ,J)*APPR!J) 

DO 62 J = 1 , N 

F(J)=C! J)*DEXP(T1!J) ) 

WR I TE ( 6, 83 ) N 

WRITE(6,R0»R01( IP) 

WPITE(6,6C ) 

WRITE(6,70)(TAU(I),F( I), EX ACT! I),I = 1,N) 

i p=ip+i - — 

R01(IP)=R01( IP-l)+0. 10+00 
IF(R01( IP). LT. 1.000) GO TO 24 
CONTINUE 
CALL EXIT 

FORMAT!//, UX, 'TIME' ,14X, ‘APPROXIMATE VALUE* , 15X, • EXACT VALUE* 



*,//) 



FORMAT! 3! 2X, 02 4. 17) ,/) 

FORMAT!//, 3X, * THE VALUE OF ROi IS , F6. 4, // ) 
FORMAT! //,3X, 'THE VALUE OF N IS *,I5.//) 
FORMAT ! //,2(2X,024. 17) ,/) 

FORMAT! //, 10 X, ‘ROOTS* ,13X, ‘WEIGHTS* , //» 
FORMAT!//, 5X,‘T *-LOG R • ,20X, / , 024. 17, / ) 
FORMAT (013.6) 

ENO 
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3. The Sphere 



The following program is listed: 

THREE-DIMENSIONAL SPHERE SOLUTION. 



The following subroutines are needed: 
INPUT subroutine 

EXPINV subroutine 
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SOLUTION OF THC- INVERSE PRORLEM IN TRANSIENT HEAT CONDUCTION 

ROUND ARV CONDITIONS ARE: 

1. INITIAL TEMPERATURE DISTRIBUTION IS ZERO. 

2. SURFACE TEMPERATURE DISTRIBUTION IS: U ( 1 . TAU I *T AU 

IV. THREE-DIMENSIONAL SPHERE SOLUTION 

DIMENSION T1!15),A(15,15>,R(15),W<15> » TAU! 1 5 ) ,C ( 1 S ) ♦ RO( 1 I I , APP 
*R l 15 ) , EXACT! 15 \ , SUM! 15 ) ,R01! 21 > 

REALMS A,R, W.CtTl, TAU, APPR, EXACT* SUM, F( 20) 

DO 100 J J = 2 » 15 
CALL INPUT (N,R,W,T1) 

CAl L EXPINVI N,R, W,A) 

- WRITE<6,106) 

DO 20 1 = 1, N 

20 WRITFI 6, 102)R( I ) ,W( I ) 

XK« 16 . OOCOOOOOOOOOOOO 
DO 4 I I = 1 , N 

4 TAIJI II) = T1 ( I I ) /XK 

PI=3. 14159265 
IR= 1 

RQ ( I R ) = 1 .OD+CO 

I P= 1 

R01I IP)=0. 00+00 
24 CONTINUE 

DO 90 !S=l,N 
SUM! I S > = 0. 00+00 
DO 6 K= 1 , 7 
S=IS+1 

6 SUM ( I S )=-SUMLI-S ) +( !-1)**k/k**3)*4 KOO+OO/ (K*#2*PT**2 + S) )*SIN(K* 

*P I *R0 1 ( IP ) ) 

APPRI IS)=( ( 1.0D+00/S**2)-! ( 1.0D+00-R01! I P ) **2) / < 6.0D+00*S ) )-(2 
*.00+00/ (PI **3# P01I IP) ) )*SUM( IS) >*R01( IP) *SINH( SORT ( S ) ) / ( S I NH { SORT! 

*S)#R01( IP) ) * XK ) 

EXACT! I S ) = T AU ( IS) 

90 CONTINUE 

00 50 I = 1 » N 

C( I >=0.00+00 

DO 50 J = 1 , N 

50 C< I )=C( I >+4< I, J)*APPR< J) 

00 62 «) = 1 » N 

62 F(J )=C( J )*DEXP(T1 ( J) ) 

WRITE(6,eC)R01! IP) 

WR I TE ( 6 , 60 ) 

WRITE! 6,70) (TA1H I ) ,F( I )* EXACT! I ) , I=1,N) 

IP=IP+1 

R01 ( IP )=ROl( IP-D+O. 10+00 
I F ( RO 1 ! IP) .LT. 1.000) GO TO 24 
100 • CONTINUE 

CALL EXIT 

60 FORMAT!//, 11X, 'TIME' ,14X, 'APPROXIMATE VALUE 15X EXACT VALUE* 

*,//> 

70 FORMAT ( 3! 2X,D?4. 17) , /) 

80 FORMAT! //^3X, ' THE VALUE- -OF— ROi-I-S *-,86. 4-,- / - / - ) 

102 FORMAT!//, 2(2X,D24. 17) ,/) 

106 FORMAT! //, ID X, 'ROOTS' , 13X, 'WEIGHTS' , //) 

108 FORMAT! //,5X, *T =-LOG R • , 20X , /, D24. 1 7 , / ) 

END — — - ■ 
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4. Subroutines 



The following subroutines are listed: 
INPUT subroutine 

EXPINV subroutine 

BESI subroutine 

& 

BESJ subroutine 



& 

Subroutine BESJ 11 and "subroutine BESI," System/360 
Scientific Subroutine Package (360A-CM-03X) Version II Programmers 
Manual . White Plains, New York: International Business Machines 

Corporation, 1967, pp. 156, 157, 159. 
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uoo I oou : i ooo’ 



r. 

c 



c 

c. 



sen Rnur INF FXP I NV{ N, R , W , A ) 

DIMENSION R ( 1 5 I , W ( l 5 ) , A (] 3 , 1 5 ) , C ( 1 7 , l 7 I , P (1 9 , j <S > , P ( 
!■' I \L*H R,W,A,C.Xl,V2,X,XI,n,P7FRn,PRCV7.PR|-Vl,/l,72 

♦ , n f \ r y 



LFf.FNIiPF COEFFICIENTS 



10 

?: 

15 

1 

AO 



#V, 



MM IIS = M4 1 - - — - 

■no 10 m=i,nplus 

C ( M , 11=1 . 01 “ no 
on ?o k=i,m>lus 
no 20 m=i,nplus 

IF(K.r,T.M) C( M t K ) =0.00+00 
CCNT INI IF 
no AO K=1,NPLUS 

nn ao v=l . nplus 

IF(K.GT.M) GO TO AO 
C(M*1,K*1)=( (2*M-1)*C(M,K+1 )-(M-l )*C(M-1,K+1 

CONTINUE 



6 COEFFICIENTS 



-no 8 f ~ I = 1 , N 

8 ( I . N ) =C ( N+ 1 » N+ 1 ) 

80 CONTINUE 

DO 90 1=1, N 

X I=R( I ) 

K = N 

50 K=K— 1 

R ( I ,K)=C(N+1»K+1 ) + X I * B ( I » K + 1 ) 

|F(k-.NE. D- G0-T0-50- 

90 CONTINUE 



LEGENDRE nF R I VAT IVES 



NV I NUS=N— 1 
DC 180 1 = 1 , N 
X = R ( I ) 

p 7 FR 0=1 .00-0 0 - 2 . on-oo*x 

PRFV 2 =P 7 ERO 
PREV 1 = 1 . 0 D 00 
nc 150 M=2 , NM I NUS 

7 2 = ?*M -1 
7 3 =y 

PRFV 3 =- ( ? 1 *PREV 1 -Z 2 *P 7 ER 0 *PREV 2 ) /Z 3 

— PRFV1 = PREV? : 

150 P R F V 2 = PR EV 3 

7 1 = N 

72=2. rn nr*x*( x-i.on on) 

130 P( I I=PREV 3 *Z 1 /Z 2 

EXPLICIT INVERSE MATRIX 



DO-220 -M=1*, N 

D F N 0 M = W < M 1 * P ( M ) 

DO 220 K=1 ,N , 

22C A(M,k)=B(M,K)/ x DENOM 

- PETIJRN 
END 




5 <? 






1 9 1 

. 73, P R F V 3 , P 



?1*C < M *x > ) / 







SliRRPUTTNF INPIIT(N,R,W.T) 

P I M|->S l ON P ( 15 > ,W< IS ), T ( 15 ) 

RFAI^-ft P.W.T 

«I-AD(5, 1 00 ) fsj 

no ir 1*1 • Mi 

10 R P A 0 < 5 * 1 0 1 ) R ( l ) 

nn li i=i,n 

11 PEAP( 5, 101 )W< n 

nn i? i=i ,n 

12 t u )=-m.nn(R( i ) ) ! 

100 FORMAT (IS) 

101 FORMAT ( 0 2 A . 17) 

RETURN 

F NO 

I 
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SUBROUTINE BESI 



PURPOSE 

COMPUTE THE I BESSEL FUNCTION FOR A GIVEN AROUMFNT AND . OR DEP 

USAGE ( ‘- 

CALL BESK x,N,BI , IER) 

DESCRIPTION OF PARAMETERS 

* -THE ARGUMENT OF THE I BESSEL FUNCTION DESIRED 
N -THE ORDER OF THE I BESSEL FUNCTION DESIRED 
B I -THF PESUI TANT I BESSEL FUNCTION 
IER-RFSULTANT ERROR CODE WHERE 
IER=0 NO ERROR 
I F R = 1 N IS NEGATIVE 

IER=2 X I S NEGATIVE 

I ER=3 UNDERFLOW, BI ,LT. l.E-69, BI SET TO C.O 
I ER=4 OVERFLOW, X ,GT. 170 WHERE X ,GT. N 

REMARKS 

N AND X MUST BE . GE. 7.ER0 
SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 

• NONE ... ... 

METHOD 

COMPUTES I BESSEL FUNCTION USING SERIES OR ASYMPTOTIC 
APPROXIMATION DEPENDING ON RANGE OF ARGUMENTS, 



-SUBROUTINE BFSI(X,N, B I , IER) — 

REAL*8 X,XX,TERM,BI ,TOL 

CHECK FOR ERRORS IN N AND X ATJO EXIT IF ANY ARE PRESENT 

I ER = 0 
B I = 1 • 0 

I F ( N ) 1 SO , IS, 1C 

10 I F < X ) 160,20,20 - - 

15 IF< X) 16C, 17,20 
17 RETURN 

DEFINE TOLERANCE 

20 T0L=1. 00-12 

IF ARGUMENT GT 12 AND GT N, USE ASYMPTOTIC FORM 

I F ( X— 12, ) AO ,40,30 
30 IF(X— FLOAT(N) ) 40,40,110 

COMPUTE FIRST TERM OF SERIES AND SET INITIAL VALUE OF THE SUM 
AO XX=X/2. 

50 TER M= 1, 0 - - - 

I F ( N ) 70,70,55 

55 DO 6C I=1,N 
F I = I 

I F(OABS( TERM) -1.0-68) 56,60,60 

56 I ER = 3 
B I =0*0 
RETURN 

60 TERM=TERMX«XX/F I - ~ — 

70 BI=TFRM 

xx= xx*xx 

COMPUTE TERMS, STOPPING WHEN ABS (TERM) LE ABSISUM OF TFRMS ) 
TIMES TOLERANCE 

DO 90 K=1,1G00 
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If (l)ABS( TH-M)-OARS(RT*Tm 1)100,100,80 

CK = K>MN«K) - - 

TFRN = 1'FRM* ( XX/FK ) 

9C HI-Rl+TERM 
C 

C RETURN MI AS ANSWER 

C 

100 RETURN 
C 

C X GT 12 -AND' X GT N,--SO-tISF— ASYMPTOT IC— APPROX IMAT ICN 

C 

11C FN=4*N*N 

IF1X-170. 0)115, in. 111 
111 I E R = A • 

RETURN 

115 XX=1./(8.*X) 

TERM-1,. 

GO 130 K-1.30 

I F ( I) A B S ( TERM)-DABS<TOL*RI ) ) 1A0, 140, 120 
.120 FK= ( ?#K- 1 ) **? 

TE«M=TFRM*XX*( fK-FN) /FLOAT ( K ) 

130 BI-BI+TERM 
C 

_ C _ SIGNIFICA NCE LOST AFTER 30 TERMS, TRY SERIFS 

GO TO AO 

1 AO PI=3. 1A1592653 

BI=BI*DEXP( X)/DSQRT(2.*PI*X) 

GO TO 100 - ' - — 

150 IFR-1 

GO TO 100 
160 I E R = 2 

- --GO-t o-l-OO 1 

END 




I 

i 
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SUBROUTINE PESJ 
PURPOSE 

COMPUTE THE J I3ESSFL FUNCTION FOR A GIVEN ARGUMFNT AND OROE 
USAGE 

CALI. BESJI X,N,BJ,D, IER ) 

DESCRIPTION OF PARAMETERS 

• X -THc ARGUMFNT OF THE J RFSSFl FUNCTION OFSIPED 
N -the ORDER OF THE J BESSEL FUNCTION DESIRED 
RJ -THE RESULTANT J BESSEL FUNCTION 
D -REQUIRED ACCURACY 
IER-RESULTANT ERROR CODE WHERE 
I E R = 0 NO ERROR 

I F R = 1 -N IS- NEGATIVE — - — 

IFR=2 X IS NFGAT I VF OR 7ERO 
I p R = 3 REQUIRED ACCURACY NOT OBTAINED 

IER = 4 RANGE OF N COMPAREO TO X NOT CORRECT ( $ FF REMARKS 
REMARKS 

N MUST RE GREATER THAN OR EQUAL TO ZERO, BUT IT MUST RE 
LESS THAN 

20 + 1 0 * X -"X * *~~2 A3 FOR— X-LESS— THAN-OR EQUAL— TO” 15 , 

9 0 + X / 2 FOR X GREATER THAN 15. 

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 
NONE 

METHOD 

RECURRENCE RELATION TECHNIQUE DESCRIBED BY H. GOIDSTFIN AND 

R-.-M.—THALERvyRECUR-RENCE— TECHNIQUES— FOR THE 'C ALCUL A T TON -OF — 

BESSFL FUNCT I0NS3,M.T. A.C. ,V, 13,PP. 102-lQfi AND I. A. STEGIIN 
AND M. ABRAMOWITZ , JDGENERA T l ON OF BESSEL FUNCTIONS ON HIGH 
SPEED COMPUTER S3 , M, T . A . C . , V . 1 1 , 1 957 , P P . 255-257 



SUBROUTINE BES J ( X , N , B J , D , I ER J 
— R E A L* 8-F M 1- , F M ,8 M K-,^d vAtF Vr1> 

BJ = .0 

IF IN) 10,20,20 

ID l ER= 1 - ' - — 

RETURN 

?0 IF( X ) 30, 30, 31 

30 I ER = 2 

RETURN : ; 

31 IFIX-15. 132,32,34 

32 NTEST=20.4lO.*X— X** 2/3 
GO TO 36 

34 NTEST=QG.+X/2. • - s t 

36 I F ( N-NTE ST ) 40 , 36 ,38 
38 IFR=4 
RETURN 

— 40- I ER -O' 

M = N+1 

6PREV=.0 
\ 

“COMPUTE STARTING VAtUE~OF-M 




50 M A = 

GC“ 

60 M A = 
70 MB = 



IF («X-5. ) 5C ,60,60 



x+6. 

TO -70 

1.4*X+6n. /X 
N4 [F IX ( X 1 / 4+ 2 



MZERO=MAXO ( MA, MB > 




SET UPPER LIMIT OF M 



MMAV=NTEST 



6Q 



r. 

c 

c 



1 1° 

1 20 
130 



140 

150 

160 
1 7 n 

ISO 

1 90 

2 00 



DO ion M=M7 FRO.MMAXj3 

St T F(M) , F ( M— 1 ) 

F M 1 = 1 .06-2 8 
F M = . 0 
4LPHA=.0 

IF ( It- (M/2 )*2) 120* 110*120 
JT = -1 

GO T0T3C " ’ 

J T = 1 
M ?. = M - 2 

CO 160 K=1*M2 

V|< _ M _ K . _ 

6VK=2.*FLCAT( MX ) *FM1/X-FM 

F V = F M 1 

FVl=8RK 

T F ( ,vk-N-T) 150Y1 40 ,150~ 

M=RMK 
J 1 =- J T 
S = I + .1 T 

AL PHA=ALPHA+RMK#S 
RMK=2.*FM1 /X-FM 
I F ( N ) 180,170*180 
8 J=RMK 

ALPHA=ALPHA+BMK 

RJ=RJ /ALPHA 

I F ( DABS( BJ-BPREV )-DABS(0*BJ ) t 200,200,190 
BPRFV=BJ 
IFR = 3 
RETURN 
ENO 
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APPENDIX B 



ROOTS OF THE SHIFTED LEGENDRE POLYNOMIALS 
AND CORRESPONDING WEIGHTS 



The roots, x^, of the shifted Legendre polynomial P^ (x) are 
computed from the roots of Legendre polynomial P^(r) using the formula 

1+r. 

1 



The coefficients (or weights) of the shifted Legendre 
* 



polynomial P^ (x) in the interval 




where A^’s are the weights of the 
interval (-1,1) and were obtained 



(0,1) are then 

Legendre polynomial P^(r) in the 
from Stroud and Secrest [20]. 
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ROOTS AND WEIGHTS 



ROOTS 



0.7 88675 134 5 948126 OD 00 
0. 211324865405186 90D 00 



0 . 88729833462074120D 00 

0 . 50000000000000000D 00 
0.11270166537925830D 00 



0.93056815579702600D 00 
0.6d999052179242760D 00 
0.33000947820757180D 00 
0.69431844202973800D-01 



0. 95 308992 2969 33150D 00 
0.769234655052S4130D 00 
0 . 50000000000000000D 00 
0.23076534494715830D 00 
0 . 46910077030668140D-01 



0.96623475710157570D 00 
0.83060469323313190D 00 
0. 6193095 930415 9840D 00 
0. 38069040695840140D 00 
0.16939530676686740D 00 
0.33765242698424040D-01 



0.97455395617137920D 00 
0. 87076 5 5 927996 9690D 00 
0.70292257568869830D 00 
0 . 50000000000000000D 00 
0.29707742431130150D 00 
0.12923440720030290D 00 . 
0.25446043828620770D-01 



WEIGHTS 



0 . 50000000000000000D 00 
0 . 50000000000000000D 00 



0.27777777777777770D 00 
0.44444444444444420D 00 
0.27777777777777770D 00 



0 . 17 3927422 56 87266 OD 00 

0.32607257743127270D 00 
0.3260725774312727CD 00 

0.17392742256872660D 00 



0 .1184634425 2809440D 00 
0.23931433524968290D 00 
0 . 2 844444444 4444430D 00 
0.23931433524968290D 00 
0 .11846 3 4425 2809440D 00 



0 . 85662246169584800D-01 
0.18038078652406920D 00 
0.23395696728634530D 00 
0.23395b96728634530D 00 
0.1803 80786 5 2406 920D 00 
0. 85662246169584800D-01 



0.64742483084434500D-01 
0. 13985 26957446 38OOD 00 

0.19091502525255910D 00 

0.20897959183673430D 00 

0.19091502525255910D 00 

0.13985269574463800D 00 

0.64742483084434500D-01 
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ROOTS AND WEIGHTS (continued) 



ROOTS 



0 . 98014492824876 800D 00 
0.89833323870661300D 00 
0. 76 2766 20495 816400D 00 
0. 59171732124782460D 00 
0.40828267875217520D 00 
0.23723379504183530D 00 
0.10166676129318d50D 00 
0 .19855071751231980D-01 



N « 8 



WEIGHTS 



0.50614268145188080D-01 
0.11119051722668690D 00 
0.15685332293894350D 00 
0. 1813418916 8918070D 00 
0.18134189168918070D 00 
0.15685332293894350D 00 
0.11119051722668690D 00 
0 . 50614268145188080D-01 



N = 9 



0 . 96408011975 3 81270D 00 
0. 918015 5 5 366 3317 40D 00 
0.80668571635029520D 00 
0.66212671170190430D 00 
0 . 50000000000000000D 00 
0. 3 3787 328829 8095 50D 00 
0.19331428364970460D 00 
0 . 619&4446336681900D-01 
0 . 15919880246 187050D-01 



0 . 40637194180787150D-01 
0.90324080347428300D-01 

0 . 1 30305 348201467 30D 00 
0 . 1561735 385 2000140D 00 
0.165H967750062960D 00 
0.15617353852000140D 00 
0.13030534820146730D 00 
0.90324080347428300D-01 
0. 406 37 194 18078715 OD-01 



N = 10 



0.986 95 3 26 42 5 85 85 40D 00 
93253168334449190D 00 
0.63970478414951200D 00 
0.71669769706462350D 00 
0. 5 744 37 16 9490615 50D 00 
0. 42556 28305091o410D 00 
0 . 2 83 3023029 3 5376 30D 00 
0 . 160295215 85 04 87 GOD 00 
0. 674683166 55507700D-01 
0 . 13046735741414190D-01 



0.98911432907302820D 00 
0.94353129988404730D 00 
0 . 86 5 0760027 8702430D 00 
0. 75954806460340560D 00 
0.6 34771577976 17220D 00 
0 . 50000000000000000D 00 
0.36522842202382730D 00 
0. 24045193 5 396 59390D 00 
0.13492399721297520D 00 
0. 5 646 8700115 95 2460D-01 
0.108856 70926 9715 20D-01 



0. 33335672154344030D-01 
0.74725674575290000D-01 
0.10954318125799080D 00 
0.1346 333596 5499780D 00 
0 . 14776 2112357 37600D 00 
0 . 147762112 35737600D 00 
0 . 1346 3335 965 499780D 00 
0.10954316125799080D 00 
0 . 74725674575290000D-01 
0 . 33335672154344030D-01 

N = 11 

0.27834283558066790D-01 
0 . 627901847 3245 1900D-01 
0.93145105463866800D-01 
0. 116 5 96 882 295 995 10D 00 
0.13140227225512.i90D 00 
0. 13646 25 433 8894990D 00 

0 . 13140227225512290 D 00 

0 . 116 5 96 882295 9 95 10D 00 
0.93145105463866800D-01 
0.62790184732451900D-01 
0. 2783428355 3086790D-01 
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ROOTS AND WEIGHTS (continued) 



ROOTS 



WEIGHTS 



N = 12 



0.99076031712335930D 00 
0.95205t>628lo523700D 00 
0.88495133709715200D 00 
0. 7936 5 897714 330S40D 00 
0.68391574949908970D 00 
0.5626167042957 3420D 00 
0. 4 373^3 2 95 74426 5 60D 00 

0.31608425050090960D 00 
0 . 206 34 102 2 85 6b 9110D 00 
0 . 115048662902847 5 OD 00 
0 . 47941371814762 c10D-01 
0.92196828766404250D-02 



0.99209152735929360D 00 
0.95879919961148860D 00 
0 . 900789045 36 b6 5 450D 00 
0 . 82 117466972016 980D 00 
0 . 72424637^ 5 1822 310D 00 
0.61522915797756730D 00 
0 . 50000000000000000D 00 
0.38477084202243230D 00 
0 . 27575 36 244 817 76 40D 00 
0 . 17882533027982970D 00 
0.99210954633345100D-01 
0.41200800388511100D-01 
0 .79084726407060300D-02 



0. 99314190434840570D 00 
0 . 96 421744183 17 86 30D 00 
0 .9136006575 3488210D 00 
0 . 84364645 240584240D 00 
0 . 75762431817907670D 00 
0.65955618446394460D 00 
0. 5 5 4027474 35 36 715 OD 00 
0.44597252564632830D 00 
0. 34044381553605520D 00 
0.24237568182092280D 00 
0. 15635354759415710D 00 
0. 86399342465117200D-01 
0. 35782558168213250D-01 
0.685 80956 515 938870D-02 



0. 23587668193255860D-01 
0 . 5 346966299765 9180D-01 
0 . 8003 916 427 1672700D-01 
0.10158371336153270D 00 
0 . 116746 26 826 9177 10D 00 
0.12457352290670110D 00 
0.12457352290670110D 00 
0 . 116746 26 82 6 9177 10D 00 
0.10158371336153270D 00 
0 . 80039164271672700D-01 
0.534696b2997659180D-01 
0.235 87668193255 860D-01 

= 13 

0.20242002382657900D-01 
0 . 46060749918864170D-01 
0.69436755109893600D-01 
0.b9072990380972600D-01 
0.10390802376844400D 00 
0.11314159013144d30D 00 
0 . 116 27 5 776 615 4 36 5 OD 00 
0.11314159013144830D 00 
0.10390802376844400D 00 
0 . 89072990380972600D-01 
0. 69436755 109893600D-01 
0.46060749918864170D-01 
0.20242002382657900D-01 

= 14 

0. 17559730165875880D-01 
O.4007904357988006OD-O1 
0.60759285 343951540D-01 
0 .786015 835 790 965 00D-01 
0.92769198738968500D-01 
0.10259923186064770D 00 
0. 10763192673157860D 00 
0.10763192673157860D 00 
0.10259923166064770D 00 
0.9276 91987 3 896 85 00D-01 
0.7 86015 835 79096 5 00D-01 
0.60759285343951540D-01 
0. 4007 904 3579880060D-01 
0. 17559730165875880D-01 
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ROOTS AND WEIGHTS (continued) 



ROUTS 



0. 99 3996 25 90102 42 20D 00 
0 . 96863669620035 2 50D 00 

0 . 92410329170521320 D 00 

0 . 86220886563008460D 00 
0 . 785 466086 3u4 2 6920D 00 
0.69707567353876130D 00 
0 .600597046 9987 1720D 00 
0.50000000000000000D 00 
0 . 39940295 300128260D 00 
0 . 30292432646121800D 00 
0.21451391369573060D 00 
0.13779113431991470D 00 
0 . 75 89670829 4786 100D-01 
0. 31363303799647090D-01 
0 . 60037409897574110D-02 



WEIGHTS 



0 . 15 3766209 9805 85 90D-01 
0 . 35183023744054020D-01 
0.5 35796102 3 35 85930D-01 
0.697853389b3076900D-01 
0.83134602908496700D-01 
0 . 93080500007780800D-01 
0 . 992157426635 55600D-01 
0.10126912096278020D 00 
0 . 9921574266 35 5 5 600D-01 
0 . 93080500007780800D-01 
0 . 63134602908496700D-01 
0. 69785 338963076900D-01 
0. 535796102335 85930D-01 
0.35183023744054020D-01 
0.1537662099d058590D-^l 
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APPENDIX C 



Legendre Coefficients 

This table gives the coefficients of the first fifteen 

& 

shifted Legendre polynomials, P^ (x) , in order of lowest to highest 
power such that 

OQ 

p / (x) -2 Tsk x< • 

K=0 

& 

where P^ (x) = P^(l-2x). The recurrence relation for these 
coefficients is 

( N -l)C Nl c = <2N-3)V lj|C - «-2>C n _ 2jK - . 

All of these coefficients are integers but due to machine conversion 
from hexadecimal to decimal, noninteger numbers may be printed in 
this table. In other words 

-0. 13999999999999990D 03 represents the integer -14. The coef- 
ficients are printed in order of lowest to highest power of x. For 
example 

P^ (x) = l-6x-6x 2 . 
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LEGENDRE COEFFICIENTS 



N = 1 

0.99999999999999900D 00 
■0 . 20000000000000000D 01 

N = 2 

0.99999999999999900D 00 
0. 59999999999999990D 01 
0. 59999999999999990D 01 

N = 3 

0.99999999999999900D 00 
0 .119 9999 9999999990D 02 
0.29999999999999980D 02 
■0 . 20000000000000000D 02 

N = 4 

0.99999999999999900D 00 
0.20000000000000000D 02 
0.89999999999999990D 02 
0.13999999999999990D 03 
0.70000000000000000D 02 



N = 5 

0.99999999999999900D 00 
•0. 299999 9999999998 OD 02 
0 . 21000000000000000D 03 
•0.55999999999999980D 03 
0.62999999999999990D 03 
■0. 25 1999999 999999 80D 03 

N = 6 

0.9 99 999999 99999900D 00 
•0.41999999999999990D 02 
0. 419999999 99999990D 03 
0.16799999999999990D 04 
0 . 314-S8999999999990D 04 
0 .277199999 99999990D 04 
0.92399999999999990D 03 



N = 7 

0.99999999999999900D 00 
-0.55999999999999980D 02 
0.75600000000000000D 03 
-0. 41999 999999999990D 04 
0.11 549999999999990D 0^ 
-0.16631999999999990D 05 
0.12011999999999990D 05 
-0.34319999999999990D 04 

N = 8 

0.99 99 99 9999999 9900D 00 
-0.72000000000000000D 02 
0.125 99999999999990D 04 
-0.923^9999999999990D 04 
0. 346 4 9 99999999 9980D 05 
-0.7 2071999 99 99999 80D 05 
0.84083999999999980D 05 
-0.51479999999999990D 05 
0.12869999999999980D 05 

N = 9 

0.999999999999999900D 00 
-0.899999999999999990D 02 
0.198000000000000000D 04 
-0.184799999999999980D 05 
0. 9008999999 99999980D 05 
-0.252251999999999990D 06 
0.420419999999999980D 06 
-0. 41183 9999999999970D 06 
0.216789999999999980D 06 
-0.486199999999999970D 05 

N = 10 

0.999999999999999900D 00 
-0.10999999 99 999 99990D 03 
0.296999999999999990D 04 
-0.343199999999999970D 05 
0.210209999999999980D 06 
-0.7 56755 999999999970D 06 
0.16 8167 999999999990D 07 
-0.233375999999999980D 07 
0. 196 91099999 99 99980D 07 
-0.923779999999999980D 06 
0.18475 5 99999 99 99990D 06 
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LEGSNDRiS COEMCIEKTS (continued) 



n - n 

0.99999999999999900D 00 
-0. 13199999999999990D 03 
0.4289 999 999 9999980D 04 
-0.b0059999999999980D 05 
0.45044999999999980D 06 
-0. 2018015 99 99999960D 07 
0.57177119999999990D 07 
-0 .105ol919999999990D 08 
0 . 1247102999 99^ 9 980D 08 
-0.92377999999999980D 07 
0.38798759999999970D 07 
-0. 705 4 3 19 9 99999 99 80D 06 

N = 12 

0.99999999999999900D 00 
-0.15599999999999990D 03 
0.60059999999999980D 04 
-0.1C009999999599980D 06 
0.90089999999999980D 06 
-0.49008959999999990D 07 
0.17153135999999990D 08 
-0. 39907295999999980D 08 
0.6235514999999996OD 08 
-0.64664599999999980D 08 
0.426786 35 999 999 990D 08 
-0. 16224935999999990D 08 
0 . 27041560000000000D 07 

N = 13 

0.99999999999999900D 00 
-0.18199999999999990D 03 
0.81O9999999999999CD 04 
-0.16015999995999980D 06 
0. 170169 999 99 999990D 07 
-0 . 11027015 99 9 999 980D 08 
0. 46558511999999980D 08 
-0.13302431999999390D 09 
0.26189162999999980D 09 
-0.35569529999999980D 09 
0.327202&7599999990D 09 
-0.19469923199999980D 09 
0.67603899999999980D 08 
-0. 10400599999999980D 08 



N = 14 

0.99999999999999900D 00 
-0 . 21000000000000000D 03 
0.10920000000000000D 05 
-0.24752000000000000D 06 
0.30630600000000000D 07 
-0.23279255999999980D 08 
0 . 116 3 96 2799999 9980D 09 
-0.39907295999999980D 09 
0. 96026 930999999990D 09 
-0.16360143799999990D 10 
0.19632172559999980D 10 
-0.16224935999999990D 10 
0 . 8768506 99999 99990D 09 
-0.2 808 16 19 9999 99 980D 09 
0.40116599999999990D 08 

N = 15 

0.99999999999999900D 00 
-0.23999999999999980D 03 
0.14279999999999980D 05 
-0. 37 127 9 9999 9999 9 80D 06 
0.52907399999999980D 07 
-0.46558511999999980D 08 
0.27159131999999990D 09 
-0.10974506399999990D 10 
0.31551705899999980D 10 
-0.6 5 44057 5 199 9 9 9970D 10 
0. 981608627999999601) 10 
-0. 1j546208399999980D 11 
0.7 9096 5 6 2999 999 960D 10 
-C.3J314267999999970D 10 
0.11633813999999983D 10 
-0.15511751999999990D 09 
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APPENDIX D 



Negatives of Logarithms of Roots of 
Shifted Legendre Polynomials 

The negative logarithms computed in these tables are 
t ± = -log x i , 

where ^ 

P N (X P = 0 ’ i = 1, 2, , N; 

N = 3, 4, ... ,15. 

It can be readily observed from these tables that the times 
or t values interlace for different N, making possible the use of 
more t values for extrapolation purposes. 
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NEGATIVES OF LOGARITHMS OF ROOTS OF SHIFTED 
LEGENDRE POLYNOMIALS 

T = -LOG(x) 



N = 2 

0.23740078615161920D 00 
0 . 15 5 4 35 86 8 307 64 310D 01 

N = 3 

0.11957401204924280D 00 
0.69314718055994510D 00 
0.21830110809448010D 01 

N = 4 

0.71959959158763400D-01 
0 . 40049171 327 565440D 00 
0.11086339030929460D 01 
0 . 2667 40966 65 219890D 01 

N = 5 

0.48046021994112870D-01 
0 . 26 2 35 92129115 56 20D 00 
0.6 9 3147 1805 59945 10D 00 
0.1466 3 5 3907 4 97 5080D 01 
0 . 30595227645462920D 01 

N = 6 

0 . 34348454520252990D-01 
0.185601297 386 3 3 380D oO 
0.4791499O102148920D 00 
0. 96576881433932520D 00 
0.17755202021841730D 01 
0.33883233221881020D 01 

N = 7 

0.2577539351C44751CD-01 
0.13838246257543370D 00 
0.35250852739216200D 00 
0.69314718055994510D 00 
0 . 12 13762486 2 357 85 OD 01 
0. 2046127 41 348787 30D 01 
0. 36711949968461940D 01 



N = 8 

0 . 20054832283440630D-01 
0.10721418964279870D 00 
0. 27080371017260680D 00 
0.5247^o25602874090D 00 
0. 8957 95 50440199 85 OD 00 
0. 14 3 87 09147125 88 30D 01 
0.22860548602906120D 01 
0. 39192958006440780D 01 

N = 9 

0.16047962730479190D-01 
0.85540945516813900D-01 
0 . 21482113 346 806 8 30D 00 
0. 41229833395649940D 00 
0. 6 9 31471805 5994 5 10D 00 
0.10850843404787260D 01 
0.16434380018108130D 01 
0.25012257285301320D 01 
0.41401866208186200D 01 

N = 10 

0.13132591978187570D-01 
0 . 6985 215 1264428800D-01 
0.174704896 357 96030D 00 
0.33310114923617220D 00 
0.55436455322540930D 00 
0.85434267909849090D 00 
0.12612407432151420D 01 
0 . 18307380647671180D 01 
0. 26 96097 17417 96480D 01 
0.4 3 39217 31125 96770D 01 
N = 11 

0.10945 35 3360154720D-01 
0 .58125740495 859540D-01 
0;i449379H41644620D 00 
0.2750316744149H30D 00 
0 . 45 4490064526 6 6760D 00 
0.69314718055994510D 00 

0.10072323073755050D 01 

0.14252350622303700D 01 
0. 200304 36 426860710D 01 
0. 28740687745821 190D 01 
0.4 5 20 307 9484 3 9 3840D 01 
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NEGATIVES OF LOGARITHMS OF ROOTS OF SHIFTED 
LEGENDRE POLYNOMIALS (continued) 



N = 12 



N - 14 



0 .926 24472047 904000D-02 
0.4912&661858812220D-01 
0.12222262181691270D 00 
0.23110140981597780D 00 
0.37992054219297830D 00 
0.57515669219094700D 00 
0.82694536135 380250D 00 
0. 11517464854216140D 01 
0.15782250281801380D 01 
0. 2 16 240008444706 20D 01 
0 . 30377764 35 14202 80D 01 
0 . 46 86 4 146 3716 346 50D 01 

N = 13 

0.79399104703031280D-02 
^.420736112037655600-01 
0.10448418267570260D 00 
0. 19701943975 657120D 00 
0 .3226236 4680743880D 00 
0.48576046598152780D 00 
0.6 9 3147 1805 5 9945 10D 00 
0.955107 33743158590D 00 
0 . 12 88247 476 8486 800D 01 
0.17213457580450050D 01 
0.23105068410166020D 01 
0.31892975959097460D Q1 
0.48398206080462550D 01 



0.68817204656915940D-02 

0.36438447740720890D-01 

0 . 90361720362083900D-01 
0 . 170021767401046 90D 00 
0.27756763865621490D 00 
0. 416 188117 9 37985 40D 00 
0 . 5905410C076837180D 00 
0.80749793054938170D 00 
0.10775051727753410D 01 
0 . 141726 635 23405600D 01 
0.18556355052547450D 01 
0.24487752135565180D 01 
0.33302947063021520D 01 
0.49823254775913270D 01 

N = 15 

0.60218359038299730D-02 
0. 31865663915715830D-01 
0 .78931426028461700D-01 
0.14825773409602940D 00 
0 . 24145 25 346 1422690D 00 
0.36086130375801570D 00 
0.50983104019704670D 00 
0 .6931471805 5 994510D 00 
0 .91778446443415 7 10D 00 

0.11942722523176470D 01 

0.15393806770033880D 01 
0.1982016 25 977 3 36 20D 01 
0.25783819644961590D 01 
0. 346 21167 382 94 8760D 01 
0.51153725057554800D 01 
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APPENDIX E 



Legendre Derivatives 

The derivatives of the shifted Legendre polynomials tabulated 
in this table are evaluated at the roots and listed in order of 
smallest to largest root. For example, 

P 5 * f (0.5000) = -3.75. 



Recalling that the recurrence relation for the derivative of 
the Legendre polynomial is 

(x 2 -l) P N (x) = N xP N (x)-P N _ 1 (x) , 
and the relationships between shifted Legendre polynomials and Legendre 
polynomials are 



and 



P N (x) = P N (1 “ 2x) ’ 

dx P N (x) = -2P N (1_2x) 



we can determine the derivative of the shifted polynomial from rhe 
formula 



-2x(x-l)^- P N (x) = N[(1-2 x)P n (x)-P n _ 1 (x)] 
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LEGENDRE DERIVATIVES 



N = 3 

-0.59999999999999560D 01 
0 .29 9999 999999999 80D 01 
-0 .5 9999999999 999990D 01 

N = 4 

0.94332756589370900D 01 
-0.37243057339H4480D 01 
0.37243057339H4550D 01 
-0.94332756589371720D 01 

N = 5 

-0 . 137 40666 961489280D 02 
0 . 4 85 177 807 2600825 0D 01 
-0.37499999999999980D 01 
0 .485177807 260083 80D 01 
-0 . 137 406 66 96 14 896 10D 02 

N = 6 

0 . 1891568687 3 5 416 6 OD 02 
-0 . 6 277070185 986 3770D 01 
0. 425 7867457 8469820D 01 
-0 . 4 257867457 8469820D 01 
0 . 6 277070185 986 4400D 01 
-0 .189 15 9 86 87 35 4217 OD 02 

N = 7 

-0. 2 4 95 6 9 9^ 2 49 8367 80D 02 
0.79712213924980360D 01 
-0.50083H9828975910D 01 
0.43749999999999980D 01 
-0 . 50083H9828975930D 01 
0.79712213924980850D 01 
-0 . 24956992249837180D 02 

N = 8 

O.318627O4825655890D 02 
-O.9923340551343249OD 01 
0.593566 167 15 355440D 01 
-0.47776387675098840D 01 
0.47 77 6 387675 098840D 01 
-0 .5 9356616715 355 900D 01 
0.99233405513433910D 01 
-0. 31862704825656630D 02 



M = 9 

-0.39632641191048770D 02 
0 . 121285105 985 90920D 02 
-0.70151138493662570D 01 
0. 5 3 499 3 85 18767 16 5 OD 01 
-0.4 9 <^18749999 999980D 01 
0.5 3499 385 1876 71720D 01 
-0.7 015 113 849 36 6 26 80D 01 
0 . 12128510598591260D 02 
-0. 3 96 326 4119105 2960D 02 

N = 10 

0.48266 5 3 972945 8070D 02 
-0.14584224552352420D 02 
0.82353816679841920D 01 
-0.60482620784350310D 01 
0 . 5 26 15 6 9 887 0016 5 80D 01 
-0.52615698870016690D 01 
0.60482620784350360D 01 
-0 . 8235 381667 9842400D 01 
0 . 145 842245 5 235 2820D 02 
-0.48266539729465570D 02 

N = 11 

-0.57764248669799700D 02 
0.17289087349651640D 02 
-0.95906545604470390D 01 
0.6 85 274800377 84 35 OD 01 
-0.57293788755183090D 01 
0 .54 1406 24 99999 99 90D 01 
-0.57293788755183200D 01 
0.685274800377o4640D 01 
-0.95906545604471760D 01 
0.17289087 34965 2200D 02 
-0.57764248669807150D 02 
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LtiGJJNDRE DERIVATIVES (continued) 



N = 12 

0.68l25674969J9b010D 02 
0.2024226 9444 5 9 3^80D 02 
0.11077665 4 369707 30D 02 
0.77531502769622o20D 01 
0.62947105594040460D 01 
•0. 57 1149 3 92612076 30D 01 
0. 57 114 9 3926 12676 5 OD 01 
•0.62947105594040800D 01 
0.77531502769623060D 01 
•0.110776654 36970840D 02 
0.20242269444595060D 02 
•0. 6812 5674989409 3 10D 02 

N = 13 

•u.793507590o 8819390D 02 
0.23443251131591050D 02 
■0 . 1269445 1159 4655 70D 02 
0.87436888274612950D 01 
•0.6941793293782O430D 01 
0 ,bll03 95609 8434 35 OD 01 
■0. 5 86 5 23437 49999S90D 01 
0.61103056098434460D 01 
•0.6 94179329 37 826740D 01 
0.87436868274013890D 01 
•0. 126 9445 115 9465 95 OD 02 
0.2344325113159191GD 02 
•0.79 35 07 5906 884 3100D 02 



N « 14 

0.9 14 3946 12 5 467 95 90D 02 
-0.26891692761957000D 02 
0.144397b9l80353500D 02 
-0.98208762794436260D 01 
0.76617260993701340D 01 
-0.6 5 88 3903 47 5 8065 5 OD Cl 
0.61321O72119930250D 01 
-0.61321072119930300D 01 
0.6 5883903475 8G6610D 01 
-0 .76o 17 26099 3701890D 01 
0.9O2J6762794437690D 01 
-0.14 4 3 976 91803 5 405 OD 02 
0. 2b 8916 9276195 92 80D 02 
-0.914394612547139 oOD 02 

N = 15 

-0 . 10439175433549440D u3 
0 . 305873 o445766o 90D 02 
-0 . 16 3127994164 5 5 820D 02 
0 . 10982495249623420D 02 
-0 . 84491321951164820D 01 
0.713285O4764050950D 01 
-0.64820514551937160D 01 
0.62841796874999980D 01 
-O.6482O514551937180D 01 
0.71328564764051460D 01 
-0.8449132l9511b5310D 01 
0.109824952496236 /OD 02 
-0 .id31279941645o730D 02 
0.30587364457665940D 02 
-0.104391754335544O0D 03 
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APPENDIX F 



Coefficients of the Polynomials 
P N (x)/(x-x i ) 

* 

The coefficients of the polynomials P^ (x) / (x-x^) are 
tabulated in order of smallest to largest root, x^, and lowest to 
highest power. For example, 

P 5 * (x) / (x-o „ 5000 ) = -i,05x o +30.4x 1 -i88.5x 2 +389.8x 3 -252x\ 

These are the coefficients of the polynomials 
* 



P N U) 



N-l 

K=0 



b .. x 
lk 



9 



* 



where x^ is a root of the shifted Legendre polynomial P^ (x) , and 
the coefficients were obtained by synthetic division 



b i,N-l ' C NN 

b i,N-2 = C N,N-l +X i b i,N-l 

b i,k-l ^N,k X i b ik 

b i,0 = C N1 +X i b il , i = 1,2, „ . . , N. 
The Legendre coefficients are defined by 



M 



COEFFICIENTS OF THE POLYNOMIALS 
P n *( x )/( x “ x i) 



N ='3 



I = 1 

-0 .11270166 5 37 9257 80D 01 
0.12254033307585160D 02 
-0.20000000000000000D 02 

I = 2 

-0 . 20000000000000000D 01 
0.20000000000000000D 02 
— 0 . 20000000000000000D 02 

I = 3 

-0. 8d72983 346 207417 OD 01 
0 . 27745966692414800D 02 
-0 . 20000000000000000D 02 



N = 4 

I = 1 

-0.107 46 123 1482126 80D 01 
0.20337454669259780D 02 
-0.74860229094208130D 02 
0 . 7 OOOOOOOOOOOOOOOOD 02 

I = 2 

-0 .14925 5 84280277 300D 01 
0 . 2762343789947&340D 02 
-0.93100663474529990D 02 
0 . 7 OOOOOOOOOOOOOOOOD 02 

I = 3 

-0.30302159969205O70D 01 
0.51422110950418290D 02 
-0.116 899 3 36 5 25 469 80D 03 
0 . 70000000000000000D 02 

I = 4 

-0.14402613260230380D 02 
0.80616996480843440D 02 
-0.13513977090579180D 03 
0 . 7 OOOOOOOOOOOOOOOOD 02 



N = 5 

I = 1 

-0.10492189413613050D 01 
0.30375707853927490D 02 
-0.18846520964827400D 03 
0.38982133941172620D 03 
-0.25199999999999980D 03 

I = 2 

-0.12999934329938890D 01 
0.37309820063988910D 02 
-0 .2 2449609 9859 34440D 03 
0.436 15 286 6 9266 8 340D 03 
-0.251999999999999SOD 03 

I = 3 

-0.20000000000000000D 01 
0.55999999999999980D 02 
-0.307 99999999999980D 03 
0.50400000000000000D 03 
-0.25199999999999980D 03 

I = 4 

-0.43334063016653720D 01 
0 .11122 377887 4647 90D 03 
-0.42803749907929230D 03 
0. 57 1847 13 307 3 315 50D 03 
-0.25199999999999980D 03 

I = 5 

-0.21317381323979400D 02 
0 . 185090693207435 20D 03 
-0.53100119141308810D 03 
0.6 1817 86605 8827 15 OD 03 
-0.25199999999999980D 03 
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COEFFICIENTS OF THE POLYNOMIALS 
P n *(x)/(x-Xi) (continued) 



N = 5 
I = 1 

■0.10349451752281860D 01 
0. 423965 85 84385 3390D 02 
■0. 39079883163O5726OD 03 
0.13342525291024040D 04 
■0.18791990844381410D 04 
0.92399999999999990D 03 
I = 2 

•0.12039421497938950D 01 
0.49116093591293450D 02 
■0 .446522767605657 30D 03 
0.14850352308906720D 04 
■0 . 20045 2126 345- 5 o50D 04 
0.9^39999999999999013 03 

I = 3 

■u . 16147012 91947 884 OD 01 
0.65210193999593780D 02 
•0.57287955811880220D 03 
0.17876688078475020D 04 
0.219975793eO29562OD 04 
0.92399999999999990D 03 

I = 4 

•0.26268064067852090D 01 
0.10342575718625110D 03 
■0 . 831579249246 3 3 320D 03 
0.22286370637292510D 04 
•0.24202420639704350D 04 
0.92399999999999990D 03 

I = 5 

■0.59033512739302930D 01 
0 . 213091197 2 416 5760D 03 
•0.12214553443508480D 04 
0.27069501770803290D 04 
•0.26154787365474120D 04 
0.92399999999999990D 03 



(N = 6) 

I = 6 

-0.29616253702314490D 02 
0. 36 6760172 137 35 010D 03 
-0. 15767 6424904777 90D 04 
0. 30574 56 19134982 80D 04 
-0. 2740800915 5 6185 40D 04 
0.923999999999999990 03 



N = 7 
I = 1 

-0 .10261104515225810D 01 
0. 5 64092826265 3 3240D 02 
-0.71785734688499980D 03 
0.35730629700534000D 04 
-0.O1852184575646220D 04 
0.86673308224198250D 04 
-0.34319999999999990D 04 

1 = 2 

-0.11484146919320910D 01 
0.62992366443543510D 02 
-0.79586014799722290D 03 
0.39093642194310200D 04 
-0.87746183861063030D 04 
0.90235324855114300D 04 
-0.34319999999999990D 04 
I = 3 

-0 . 14226 31787035 9660D 01 
0. 77643498872 5 338 IOD 02 
-0. 96505152144649300D 03 
0.46021405350141410D 04 
-0. 98842457267481930D 04 
0.95995697202363800D 04 
-0.34319999999999990D 04 
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CO aPiU Cl MINI’S OP 
P n *(x)/(x-x. ) 



(N = 7) 

I = 4 

■0.20000000000000000D 01 
0. 107 9 9999 99999 9990D 02 
■0. 1295 999999999 9990D 03 
0.58080000000000000D 04 
•0.1148 3 999999 999 9 80D 05 
0.10295 999 9999 9 99 90D 05 
•0.34319999999999990D 04 

I = 5 

•0. 3 365105 58665 4900D 01 
0.177172244788888801) 03 
•0 . 194 840707 4 529 3940D 04 
0. 7579145 16 96147 850D 04 
•0.13366397125566240D 05 
0 . 109 924 3027976 3600D 05 
■0.34319999999999990D 04 

I = 6 

•0.77378774094586310D 01 
0.37344638812587250D 03 
•0.29601529512276080D 04 
0 . 95 9 37 84 4 6987 97 980D 04 
•0.1513695595O549090D 05 
0. 1156 84675 144oo540D 05 
•0.34319999999999990D 04 

I = 7 

392988 3 98013 854 20D 02 
0.6 56 3 36 2191426 2 670D 03 
G.39166709579142760D 04 
0 . 1113 450263 6006 810D 05 
0 . 16 32 856 4 34 546 5 460D 05 
0.11924669177580150D 05 
•0.34319.999999999990D 04 



THE POLYNOMIALS 
(continued) 



M = 8 

1 = 1 

-0.10202572815328950D 01 
0. 7 2417 5 99 349605 4 80D 02 
-0. 12116 395 9167 9 8290D 04 
0. 619099 3165 3 3167 20D 04 
-0. 26 994994 385 108750D 05 
0. 45 99014 3208138260D 05 
-0. 3 8865 5 34 77 34 383 50D 05 
0.12o69999999999980D 05 

I = 2 

-0.11131726590019590D 01 
0.78909276079304370D 02 
-0.1 3147 57 89944 15 500D 04 
0. 882 216 2 822 30317 20D 04 
-0. 267 50o4 2167 3346 50D 05 
0. 482^39284 33 314890D 05 
-0. 3 99184 51217 843 300D 05 
0. 1266 9999999 9999 80D 05 

1=3 

-0.1311017705o897280D 01 
0 . 92674507 384 851140D 02 
-0 .15303843891165150D 04 
0.10107442569908190D 05 
-0. 32175727 33 3695730D 05 
0.52304719856449930L 05 
-0 . 41663198942188390D 05 
0.12669999999999980D 05 

I = 4 

-0 . 16 8999615 87 9 37 5 20D 01 
0.118823636416346oOD 03 
-w . 19265836709615 300D 04 
0.123562655H410070D 05 
-0 . 3 7 676 3 25 65086 3 890D 05 
0.58128557529128660D 05 
-0 . 4 38645 98075 540460D 05 
0.12 66 9999 99 9 9999 80D 05 
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COEFFICIENTS OF THE POLx KOIV.IALS 



P n *(x)/(x-x i ) 



(N = 8) 

I = 5 

-0.2449263430431910UD 01 
0.170349417o6849940D 03 
-0 . 266 886 3 116 26487 20D 04 
0 . 1609457 66 8843 14 80D 05 
-0.454474908616721000 05 
0 . 65 210969075 865700D 05 
-0 . 46 2 25 40192445 9470D 05 
0.128699999999999800 05 

I = 6 

-0.42152510346328980D 01 
0.2O572973320859910D 03 
-0 . 41067 9 375 0104 5 800D 04 
0.21O37752955856010D 05 
-0.548490878157275000 05 
0.7259552d2C3319190D 05 
-0.48426801057811540D 05 
0.12869999999999980D 05 

I = 7 

-0 . 98360564168348000D 01 
0.6 1144605 6 2816 07 80D 03 
-0. 63791935 -69591110D 04 
0.26139054049247060D 05 
. 64042031711590040D 05 
0.78983221126254860D 05 
-0 . 501715 4878215 6 66 OD 05 
0.12669999999999980D 05 

I = 8 

-0.503649653H09^020D 02 
0.106964777161117300 04 
-0. 857 978407 346 99980D 04 
0.33251752237512180D 05 
-0.7 04 ^270005 4007100D 05 
0. 83*j66934567508040D 05 
-0 . 5 1224465 2 26 5616 10D 05 
0. 1286 9999 999999980D 05 



( continued) 



N = 9 

I = 1 

-0.101617742287028100 01 
0.90423351504540900D 02 
-0.192014512900450900 04 
0 . 16 82775 06 4609 3990D 05 
-0.74447443742928090D 05 
0.160680975753836800 06 
-0 . 24361738382250640D 06 
0.1709440245775 6 9 30D 06 
-0 . 486 199999 9999 9 97 OD 05 

I = 2 

-0. 1089 3061626 26 94bOD 01 
0.968509667211815600 02 
-0.205132584711027900 04 
c. 176956559966609200 05 
-0.78641525966798870D 05 
0 . 1891149592 d2305 40D 06 
-0.25196200632405130D 06 
0. 17415 60837 80889 30D 06 
-0.48619999999999970D 05 

I = 3 

-0.12396401469971320) 01 
0.11003090553602840D 03 
-0 . 23180887631484420D 04 
0 . 2003495 402 146 6 020D 05 
-0.65643047494983990D 05 
0.205047578198581400 Oo 
-0.26698430062184840D 06 
0.179568940471048300 06 
-0.486199999999999700 05 

I = 4 

-0.15102649384090920D 01 
0.133644683861402300 03 
-0.27885226249108700D 04 
0.23o98b01940943240D 05 
-0.10026992852834280D 06 
0.22953623345146140D 06 
-0.2 882 68877 6045 5 7 10D 06 
0.16659739927705310D 06 
-0.486 199999999 99970D 05 
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COEFFICIENTS OF THE POLYNOMIALS 
P n *( X )/(x- Xi ) (continued) 



(N = 9) 

I = 5 

-0 . 20000000000000000D 01 
0. 17599999999999990D 03 
-0.36079999999999980D 04 
0.29743999999999990D 05 
-0 .120691999 999999 80D 06 
0.263H999999999980D 06 
-0. 31459999999999980D 06 
0. 19447999999999970D 06 
-0.48619999999999970D 05 

I ='6 

-0. 295 96 894298 360320D 01 
0. 257o 12287 16 4207 20D 03 
-0.50977327077604850D 04 
0. 39607 35 505 2089010D 05 
-0.14941295064252480D 06 
0.30437164735776020D 06 
-0.3434670826651C310D 06 
0. 202 36 260072 294640D 06 
-0.48619999999999970D 05 

I = 7 

-0.51729234959793730D 01 
0.438803977142&9750D 03 
-0.7972489H85476270D 04 
0.54354549922924990D 05 
-0.18485674934310280D 06 
0. 34863047564048830D 06 
-0.37136171732450760D 06 
0.20939105952895120D 06 
-0.48619999999999970D 05 

I = 8 

-0 . 121974 35 5 464 125 70D 02 
0.94899176526824240D 03 
-0 . 12 575 65c 49096 1220D 05 
0.7 2017 849 395 379870D 05 
-0.22043389218491320D 06 
0.38809931928332460D 06 
-0. 39422941985782490D 06 
0.21480391621911020D 06 
-0.48619999999999970D 05 



(N = 9) 

I = 9 

-0.6 28145 4285 6 8479 30D 02 
0. 17076 4206 28015 200D 04 
-0.17108039318556490D 05 
0.86179083022441700D 05 
-0.24566246209640380D 06 
0.41391861105224000D 06 
-0.40836921177951900D 06 
0.21801597542243020D 06 
-0.48619999999999970D 05 



(N = 10) 

I = 1 

-0.10132192030187380D 01 

0.11042749919759750D 03 
-0.28973737707332620D 04 
0.31O36008310223160D 05 
-0.18072992729171640D 06 
0 . 5 836 4067 8407404 20D 06 
-0 .11125 5 4526 4977240D 07 
0. 12 37 3488 36795 895 OD 07 
-0 .7414 3446 27086 4040D 06 
0. 18475599999999990D 06 

I = 2 

-O.1072349623863344OD 01 
0 . 116 8085 249 -3755 10D 03 
-0. 30596188054902030D 04 
0.33522058019943660D 05 
-0.18947124814716320D 06 
0. 6083275 9034876 93uD 06 
-0.11510090529060490D 07 
0. 12683225333985980D 07 
-0.7 5 14 89176 3 12004 80D 06 
0.18475599999999990D 06 
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COEFFICIENTS OP TflE POLYNOMIALS 

* 

P n ( x )/( x “ x i) (continued) 



( N = 1 0) 

I = 3 

0. 11908947 27 30826501) 01 
0.1295 8018976025 95 OD 03 
0 . 3 3826 40975 5 6 2185 OD 04 
0 . 368431 3774 1286 600D 05 
0.20646168216643710D 06 
0.65 5 342601615 54700D 06 
u . 12 222 5 97 96 2496 65 OD 07 
0.13236797321288390D 07 
•0.76 86 39 5 02 89 967 240D 06 
0.18475599999999990D 06 

I = 4 

•0. 1395 26642 36 969410D 01 
0.15153469oo2067490D 03 
•0.39325717310421160D 04 
0. 42 399226889 406 180D 05 
•0.23414442909178530D 06 
0. 72919387497500 370D 06 
•0.13289928639733150D 07 
0.140193995 34028090D 07 
•0.7 9136 5 8002 8112 810D 06 
0.18475 5 99999 999990D 06 

I = 5 

■0.17408344256148520D 01 
0. 18846 1282 3197 8200D 03 
•0. 4842198 35 5 906 85 90D 04 
0.51315971893362J30D 05 
•0.27660819415199490D 06 
0. 83585782980166610D 06 
■0.14724363518260360D 07 
0. 149942185 8 33 385 40D 07 
•0.81764928d31355450D 06 
0. 18475599999999990D 06 



(N =10) 



I = 6 

-0. 2 3 4982926091 55 160D 01 
0.25295952325131730D 03 
-0.63845812696981660D 04 
0.65643464907H4140D 05 
-0. 3397066772 J184200D 06 
0.97999471029732480D 06 
-0 . 16488406397 314120D 07 
0 . 160944 356 782 5 4130D 07 
-0.84515471366644480D 06 
0.18475599999999990D 06 

I = 7 

-0.3 5 2979 834487 39 9 30D 01 
0. 37561834158055780D 03 
-0.91569381241888310D 04 
0.68820534161174900D 05 
-0.4 284803 35 60290120D 06 
0. 115 87 46 896 85 096 20D 07 
-0.18458484019747700D 07 
0.17 222295 5 115 37820D 07 
-0. 8714 381997 1887 HOD 06 
0.18475 5 99999999990D 06 

I = 8 

-0 .6 236489 36909 479 10D 01 
0.64731508099210330D 03 
-0. 14490045174987220D 05 
0.12370896236546900D 06 
-0.53963580369c67550D 06 
0 . 135 45 02036 441 5 5 10D 07 
-0.20410962474613890D 07 
0. 1825 77 97089 3145 60D 07 
-0 . 89416449710O3271OD 06 
0. 16475 599999999990D 06 



82 



COEFFICIENTS OF THE POLYNOMIALS 
P n (x)/(x-x^) (continued) 



(N = 10) 

,1 = 9 

0. 14S21771900816640D 02 
0.14107099868099880D 04 
0. 231114409027236 30D 05 
0.16613070627655260D 06 
0.65333323711654740D 06 
0.15329086008107720D 07 
0.22050557441480030D 07 
0.1907625122y025560D 07 
•0.9H31482368799460D 06 
0.18475 5 999 99999990D 06 

I = 10 

•0. 766475 24700426040D 02 
0.25 56 384674 344 345 OD 04 
•0 . 3 17025 908896 6705 OD 05 
0.2006 17 9 294 35 367 10D 06 
■0.735 2084655 3093450D 06 
0.16515651804509890D 07 
•0.23082263752316030D 07 
0. 195 7089135 12676 80D 07 
•0.92136953729135910D 06 
0.18475599999999990D 06 



N = 11 ' 

I = 1 

0.10110054727960090D 01 
0.13243059035446710D 03 
O.42033254270432940D 04 
0.56471403690313600D 05 
■0.39831451706792320D Oo 
0.16.375270636613190D 07 
0.41250893009936650D 07 
0.64470107363448310D 07 
0.60903164443090320D 07 
0.31821231006133540D 07 
•0.705 4319999 9 99 9980D 06 



IN fJD 
I = .2 

-0 . 10596482530719400D 01 
0.13877669109972230D 03 
-0. 4 39966677 2485 89 90D 04 
0.58991506942434560D 05 
-0.41468660005838970D 06 
0.16990738941449220D 07 
-0.42591465766413190D 07 
0.66163925077269420D 07 
-0.62050273191705860D 07 
0.32142788280601930D 07 
-0.705 4 31999999 9 99 80D 06 
I =_3 

-0.115 5 96779571 35 300D 01 
0. 15123148747941770D 03 
-0.47642599947134420D 04 
0. 63896975 326 10325 OD 05 
-0.44664264782901690D 06 
0.18162255652787900D 07 
-0. 45099926736514970D 07 
0.69264750230549010D 07 
-0 .6 409 326994486 1990D 07 
0.32696237052019400D 07 
-0.705 4 31999 99 999980D 06 

I = 4 

-0.13165723758365340D 01 
0.17205419079357020D 03 
-0. 54215736977179900D 04 
0.71935442730423050D 05 
-0. 49834180996460860D 06 
0.20007610589184740D 07 
-0. 489 36 34 9314803 5 OOD 07 
0.73837131971996860D 07 
-0.66978207698503270D 07 
0.33440664896906870D 07 
-0.70543199999999980D 06 
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COEFFICIENTS OF THE POLYNOMIALS 

41 

P n ( x)/(x-x i ) (continued) 



(N = n ) 

T «= 5 

0.15753698412222740D 01 
0.20546702890291260D 03 
•0. 64 3465 005 808815 40D 04 
0.84479759022740410D 05 
■0.576 5 384804 15 366 10D 06 
0 .22708602111 336490D 07 
■0. 5 430066 35 5 2 9565 30D 07 
0.79900452708907000D 07 
■0 .7 05 920820115 46 810D 07 
0. 34320878162051090D 07 
•0.70543199999999980D 06 

I = 6 

•0 . 20000000000000000D 01 
0.25999999999999980D 03 
•0.8O59999999999998OD 04 
0.10399999999999990D 06 
■0.69289999999999980D 06 
0.26502319999999990D 07 
•0.61349599999999980D 07 
0.87339199999999990D 07 
•0.74742199999999980D 07 
0.35271599999999980D 07 
•0.70543199999999980D 06 

.1 = 7 

•0.273801254146d0460D 01 
0.35392094279590580D 03 
■0 .1077703 382 2814850D 05 
0.13493737947363220D 06 
•0.86387751199106730D 06 
0. 31600456 5 47 5365 80D 07 
•0.70029225301624410D 07 
0.95 8029895 496 27140D 07 
•0.79148578553086700D 07 
0.36222321837948840D 07 
•0.70543199999999980D 06 



(N ■ 11j 

• 1*8 

-0.41588353129710050D 01 
0. 5 3167035 015 1982 20D 03 
-0 .15 630274065 5 80460D 05 
0.18477591316176 540D 06 
-0.11048947740846370D 07 
0.37975207993617690D 07 
-0. 7 9 85 75 8972 873 82 90D 07 
0.10464299332738090D 08 
-0.83456623626341160D 07 
0.37102535103093070D 07 
-0.705 43199999999980D 06 

1 = 9 

-0.7 4115 80005 456 6 390D 01 
0.92339704254301420D 03 
-0. 24951847165799960D 05 
0.26020o88357447880D 06 
-0 . 14100020778751850D 07 
0.45063438282596450D 07 
-0.89781521209175620D 07 
0.11293527545564720D 08 
-0. 87271536476687060D 07 
0.37846962947980530D 07 
-0.70543199999999980D 06 

I = 10 

-0. 177089 25 439165 380D 02 
0.^U239721177599240D 04 
-0.40128918809659610L 05 
0.35295803072169190D 06 
-0.17264780148669660D 07 
0 . 5 16 2 824441405 3760D 07 
-0.98264624022726930D 07 
0.11961628235470690D 08 
-0. 9020957 8666 288160D 07 
0. 38400411719398010D 07 
-0.70543199999999980D 06 
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COEFFICIENTS OF THE POLYNOMIALS 
P n (x)/(x-x ± ) (continued) 



( N - 1 1) 

I = 11 

■0.91863882962169210D 02 
0.36870595581192030D 04 
■0.55388450186096020D 05 
0.429146705 35 6415 70D 06 
■0.19570033658468080D 07 
0. 5 6048 5746 3062 36 20D 07 
•0. 103672 5 4135 710960D 08 
0.12370929196046650D 08 
■0.919 56485 3 87888340D 07 
0.387 2196 899 3866 440D 07 
•0.70543199999999980D 06 



N * 12 



I - 1 

•0.1009 305 4785 0814 20D 01 
0.15 64 3 2 95 6774406 90D 03 
■0.59040000514031970D 04 
0.95072538604810130D 05 
■0. 81332606983 4 36670D 06 
0. 4125 6067 15 7 5 86 6 5 OD 07 
•0.13148756650783650D 08 
0. 27007540306116740D 08 
•0.35676536042330500D 08 
0.29257811703237760D 08 
■0.1354571146C768950D 08 
0 . 27041560000000000D 07 

1 = 2 

■0.10503554820636720D 01 
0.16275220866733050D 03 
■0.61374873 5 46 0835 30D 04 
0.98694040328690540D 05 
•0.84260142802385u90D 06 
0 . 426265 195 42 3085 40D 07 
0.13539590592588070D 08 
0. 27695 46 3941829490D 08 
■0.36404991281091360D 08 
0.2 96 826 34957 866470D 08 
■0 . 1365042O94824111OD 08 
0 . 270415 60000000000D 07 



(N - 12 ) 

1 = 3 

-0. 11300056386074940D 01 
0.17 5 0039667 37 20190D 03 
-0.65890583909278350D 04 
0.10566789120383710D 06 
-0. 89861676621079580D 06 
0 . 45225 980977865030L 07 
-0.14272579036542980D 08 
0.28967374689262430D 08 
-0.37728374331019390D 08 
0. 30438086864005020D 08 
-0.13831889532080700D 08 
0.27041560000000000D 07 

I = 4 

-0.12599870085909810D 01 
0.19497040599021400D 03 
-0.7 321 821791679 3060D 04 
0.11689929917036390D 06 
-0.9878 3069732 489 310D 06 
0.49304114429093590D 07 
-0.15 4004741445 56 430D 08 
0. 3087 8277145 7499 80D 08 
-G.39660450849491590D 08 
0. 315 0490307626 5 990D 08 
-0.14078758315004040D 08 
0.27041560000000000D 07 

1=5 

-0. 146 2lb 84041825 5 10D 01 
0 . 2 ^596033460115 980D 03 
-0. 845 1391373 326 65 10D 04 
0 . 134005699815 80050D 06 
-0.11213286150317220D 07 
0.55263640115562260D 07 
-0.17000298643450430D 08 
0. 3 349388776 8086 840D 08 
-0. 4220002 57 3 8625 690D 08 
0 . 32 8469906 97067070D 08 
-0.14375521122497510D 08 
0. 270415 bOOOOOOOOOOD 07 
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COEFFICIENTS OF THE POLYNOMIALS 
P n (x)/(x-x i ) (continued) 



(N = 12) 



1 = 6 

0 . 17774090110084440D 01 
0.274116622955u0000D 03 
0 .10187901165 5 49HOD 05 
0.15 98105746 849 3 100D 06 
0.13172190226655770D 07 
0.63696597527710040D 07 
0. 19 166 6 47 85 396 66 90D 08 
0. 36 8646 14 91 304 336 OD 08 
0.453071067640572000 08 
0.344061829119526100 08 
■0.14703532663486610D 08 
0. 270415 50000000000D 07 

1=7 

■0. 2 286 32416 85 604 96 OD 01 
O.351439292097O5O4OD 03 
0.12928158809268100D 05 
0.199303087335325400 06 
0. 16040779760371220D 07 
0.75375947230744230D 07 
0.21984244415561060D 08 
0.40977 906 0124 83700D 08 
0.48875309586617870D 08 
0.36099436277086250D 08 
■0. 15 042183 3 36 51335 OD 08 
C.27041560000000000D 07 

1 = 8 

0.31637134669490470D 01 
0.48353021794298810D 03 
0.17471512020309010D 05 
0.26 141286017 4 92700D 06 

0 . 20231540762048540D 07 

0 . 9104 3 5087 87125 C80D 07 
0.25464049880790580D 08 
0.45694292253791910D 08 
0.52710180022588890D 08 
0. 37 82035 947 2091660D 08 
0.15 370194 877502440D 08 
0.27041560000000000D 07 



(N = 12 ) 

1 = 9 

-0. 4 846 34604479026 40D 01 
0.73254291300179080D 03 
-0.25556997895957640D 05 
0.3ol26ll8341389720D 06 
-0.26152764443786540D 07 
0. 11076 903 29 3286 250D 08 
-0.29447526345421890D 08 
0.50691663294906840D 08 
-0.56525292661719900D 08 
0. 39445 899926 225 400D 08 
-0.156669576 84995 930D 08 
0 . 27 041560000000000D 07 

I = 10 

-0.86919741157220160D 01 
0.12803975480242770D 04 
-0 . 4 1074814 19 37 6626 OD 05 
0.51304538720347610D 06 
-0.3371222255o908350D 07 
0. 13^95 8845 96 2 5026 OD 08 
-0.33527129359225730D OC 
0. 5 5 456243295 603750D 08 
-0. 59965118501394010D 08 
0.40847771543197920D 06 
-0.15 9 13 826467 919280D 08 
0.27041560000000000D 07 

I = 11 

-0. 20858810712881270D 02 
0.28188844868536790D 04 
-0.66479439208806720D 05 
0.70128491360441830D 06 
-0.416 37333025608900D 07 
0.15376337170479350D 08 
-0.37061910459839990D 08 
0. 5 935 135 83 873666 90D 08 
-0.62655520668857870D 08 
0.41907005475455270D 08 
-0.16095295051758860D 08 
0. 27041560000000000D 07 
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COEFFICIENTS OF THE POLYNOMIALS 
P n (xJ/Cx-Xj^) (continued) 



( N = 1 2) 

I « 12 

0 .10846 360047086 500D 03 
0.5 15 5 96 9046 35 14 85 OD 04 
0.92197417744400200D 05 
O.85714252445951830D 06 
0. 474609 3 3446 36 4170D 07 
0 . 1679045 3 36 3184 5 OOD 08 
C.3 93 37 8646 1727 2 180D 08 
0.6 176 2577 99175 76 80D 08 
0.6 4 27 249 35 5 2206 250D 08 
0.42529277095548140D 06 
•0.16200004539231030D 08 
0.27041560000000000D 07 



N « 13 
1 = 1 

0.10079715286689510D 01 
0.1624348091o835120D 03 
0.607139761o0663870L 04 
0. 153 30097 897 99 9020D 06 
•0. 15 607 42107 2745 86 OD 07 
0.95417344384769900D 07 
•0.373H857364665380D 08 
0. 96475 435 981292820D 08 
•0.16 67 348117 15 41170D 09 
0. 1904 26 470819127 40D 09 
•0.13786672036695810D 09 
0.57285552860546860D 08 
•0.104005 99999999980D 08 



(N - 13 ) 

1 = 2 

-0.10429712616983220D 01 
0.18873297851273380D 03 
-0.83450914693393850D 04 
0.15 83385 8496 35 4280D 06 
-0 . 16096 815 846 96 6000D 07 
0.98220090495688290D 07 
-0.38315116413652590D 08 
0. 98778976478833180D 08 
-0.17012180818179740D 09 
0.19350609793310380D 09 
-0.13944189578075440D 09 
0.57631813044520720D 08 
-0.104005 99999999 980D 08 

1 = 3 

-0 . 11101378477288740D 01 
0.20081267981965990D 03 
-0. 886 909910957 947 20D 04 
0.16795375306639010D 06 
-0.17026697369630140D 07 
0.10351309566871580D 08 
-0. 401949 85 2957 977 50D 08 
0.10305 33 5 65 896 2 990D 09 
-0. 17633237685047210D 09 
0.19907316154864720D 09 
-0.14224164371270630D 09 
0.58235153454759510D 08 
-0.104005 99999999980D 08 

1 = 4 

-0. 12177 67 7100408040D 01 
0.22015076567274520D 03 
-0.97054250797130950D 04 
C.18321872369925510D 06 
-0. 184 915747196 41470D 07 
0. 11176 4997 95 303420D 08 
-0.43067072104590820D 06 
0. 1095 226767 3582 370D 09 
-0.18554999183803190D 09 
0. 207 1487 5 218926 85 OD 09 
-0. 146 19 8035 85 95 836 OD 09 
0.59063190730108370D 08 
-0.104005 99999999980D 06 
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COEFFICIENTS OF THE POLYNOMIALS 
P n (x)/(x-x^) (continued) 



(N * 1 3) 

1 = 5 

-0.13807456076844240D 01 
0.24938924169703470D 03 
- 0.1096 396 3406 2665 10D 05 
0.20600177182381960D 06 
-0.2065178755O96915OD 07 
0. 12 374017 3894 4 2320D 08 
-0 .472000906 9579045 OD 08 
0.11850142742212430D 09 
-0. 19798539202253470D 09 
0.21770203249501490D 09 
-0.15119280842328430D 09 
0.60071303146785120D 08 
-0.10400599999999980D 08 

1 = 6 

-0.16254106083538880D 01 
0. 29 318277109067 5 40D 03 
-0.12835570496802190D 05 
0.239 4626905 9726430D 06 
-0.23767360347638980D 07 
0. 1406025 682149 39 90D 08 
-0.52823008723021120D 08 
0. 1 303600621606 15 30D 09 
-0.21379280571123440D 09 
0.23058480315709960D 09 
-0. 1570440405 6613840D 09 
0.61205147619538480D 08 
-0.10400599999999980D 08 

1 = 7 

-0 . 20000000000000000D 01 
0.35999999999999990D 03 
-0.15 65 999999 9 999980D 05 
0.28899999999999990D 06 
-0.2825 399999999998OD 07 
0.16403231999999980D 08 
-O.6O310559999999990D 08 
0.14542751999999990D 09 
-0.23292d21999999980D 09 
0.24545415999999970D 09 
-0.16349743199999980D 09 
0. 624035 99999999980D 08 
-0.10400599999999980D 08 



(N = 13) 



1 = 8 

-0.25989495325629240D 01 
0.46625427624523020D 03 
-0.2007 3625338025120D 05 
0.36407741793960490D 06 
-0.34764135843287480D 07 
0. 196 2 36 346 1842 316 OD 08 
-0.70002386979225560D 08 
0.16379082336259800D 09 
-0.25495904554972060D 09 
0.2617O448342966820D 09 
-0.17022701675121470D 09 
0.6 360205 2 380461440D 08 
-0. 104005 99999999980D 08 

1 = 9 

-0. 3626 425 5 8872 2326 OD 01 
0. 646 85 84945 9456770D 03 
-0.27354641374456O7OD 05 
0.48160875083736360D 06 
-0.442 4570126 5 2 34570D 07 
0 . 2394327866 3642210D 08 
-0. 82012460865595080D 08 
0.18499071129263050D 09 
-0.27887545939564370D 09 
0. 27 8436 3786 646 4460D 09 
-0.17684807380864750D 09 
0.64735896853214790D 08 
-0.104005 999 9 99 99980D 08 

I = 10 

-0. 5 5 920489476 264 360D 01 
0 . 98648189703509230D 03 
-0. 40282425 826 881880D 05 
0. 67036126 249861210D 06 
-0.57672967016894390D 07 
0 . 2 9412607766 919960D 08 
-0. 95 8807 35722672 390D 08 
0 . 20770874136904660D 09 
-0.30299336534800010D 09 
0.29448811625060740D 09 
-0.182 9425 37 82839130D 09 
0. 65744009269691 540D 08 
-0. 104005 99999999980D 08 
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COEFFICIENTS OF THE POLYNOMIALS 
? n (x)/(x-x i ) (continued) 



( N = 13 ) 

I * 11 

O. 100795 32080864670D 02 
0.17 32877 8717 4 816 80D 04 
■0. 6 508476 964177 47 20D 05 
0.9p831383449131930D 06 
•0.74929847037156350D 07 
0 . 35 6 2 13818256775 50D 08 
0 . 11024 115 446 467 420D 09 
0.22964364791696370D 09 
•0. 32504456994910890D 09 
0.30854183556664090D 09 
•0.18809455571035100D 09 
0.66572046545240380D 08 
•0 . 104005 99999999980D 08 

I = 12 

•0.24271373142518990D 02 
0 . 38282903577150800D 04 
•0. 105 8646 822 672 2040E 06 
0.13178219165 5 39170D 07 
•0.93172482038267870D 07 
0.414984 1216 79825 10D 08 
•0 . 12281557116129330D 09 
C.24778035238250510D 09 
•0.34250008457189330D 09 
0.31929514242560010D 09 
•0.191931552291026100 09 
0.67175 3 86 95 547 9200D 08 
•0.104005 99999999980D 08 

I = 13 

•0. 12644666618088360D 03 

0.70245 3385 66 609130D 04 
•0.14736930837191630D 06 
0.16 17 340314 5 188670D 07 
•0.10667000989156780D 08 
0.45 5 22697896 197000D 08 
•0.13097524020902070D 09 
0.25909930830793510D 09 
•0.35 307976 886 615020D 09 
0.32566732552037570D 09 
•0.19416523890094200D 09 
0.67521647139453030D 08 
•0.104005 99999 9999 30D 08 



N = 14 



1*1 

-0. 1006 905 478 2458 36OD 01 
0.210436 2867045 3 990D 03 
-0.10783518111968030D 05 
0. 3837125 47 5 3723 30D 06 
-0. 2 844194503 2009650D 07 
0.20576174872216570D 08 
-0. 96 4 817 8644787 95 90D 08 
0.304680702956188900 09 
-0. 6601157 4 395 699110D 09 
0. 982b 3765 909997460D 09 
-0.98735094411646760D 09 
0.63952860422322550D 09 
-0. <^4097472348001660D 09 
0.40116 5 99999999990D 08 

1 = 2 

-0. 1037 11046 84 22 8440D 01 
0. 2 167175 99647022 30D 03 
-0.11100486193258o90D 05 
0.24519315203176540D 06 
-0. 292 24 391985 846 560D 07 
0.21112267750251590D 08 
-0. 98820046 30483 6740D 08 
0.3H39543910837490D 09 
-0.672953&82330646200 09 
0.998600121101071300 09 
-0.10002071037698320D 10 
0.6 45 37 983781746 240D 08 
-0.24213507457301090D 09 
0.40116599999999990D 08 
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COEFFICIENTS OF THE POLYNOMIALS 
P n *(x)/(x-x ± ) (continued) 



(N - 14) 



■ 1 = 3 

0.10945701042210040D 01 
0.22 86 61645 7 34205 90D 03 
0. 117024197236 39010D 05 
0.258118881955H410D 06 
•0. 307020479 34 305 350D 07 
0.22120224016802250D 08 
■0.10319175583517800D 09 
0. 32386273H8846210D 09 
•0.69659163832994400D 09 
0.10282640822576100D 10 
■0.10233718266634350D 10 
0.6 5 57808035 6 5 904 50D 09 
•0.24416564786193610D 09 
0. 40116 599999999990D 08 

I = 4 

■0. 1185 3306 3847006 60D 01 
0.24751442830826710D 03 
■0 . 126504242876 34700D 05 
0.27839810745672320D 06 
•0.33007450983788320D 07 
0.23681141365139410D 08 
■0. 109898095 7846 2 980D 09 
0.34276783051801440D 09 
0.73194343166031770D 09 
0.10716230071987210D 10 
■0.10568339927924010D 10 
0. 670493 67136484 310D 09 
■0.24697197272741550D 09 
0.40116 5 99999999990D 08 



( N « 1 4) 

1*5 

-0.131991539o2437980D 01 
0.27544005597036710D 03 
-0.140499185 26 3 92 3 90D 05 
0.30816075444191710D 06 
-0. 36 362339215 55 4540D 07 
0.25927127214786090D 08 
-0.11941162739159840D 09 
0.36912929785643290D 09 
-0.7 802 5 48016 995 4 270D 09 
0.11295302404722740D 10 
-0.11003963250961380D 10 
0.68912423001244720D 09 
-0.25042288827733710D 09 
0. 401165 99999999990D 08 

1*6 

-0 . 15 16 17 10607 407 5 80S 01 
0.316097148128032100 03 
-0 .16 077 3 306 37 8596 10D 05 
0.350906677571745700 06 
-0.41120883805108430D 07 
0.290607048663598800 08 
-0.13241567161503020D 09 
0.40429806385895040D 09 
-0.84294751415744180D 09 
0.1202425031443H00D 10 
-0.11534911543210250D 10 
0.711087935685351400 09 
-0. 25 43570483703 3 370D 09 
0. 40116 5 999999999 90D 08 
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COEFFICIENTS OF THE POLYNOMIALS 
P n *(x)/(x-x i ) (continued) 



( N = 1 4) 

1 = 7 

■0. 1804 96464-0 35 377200 01 
0. 37 57 8467674825 4 20D 08 
•0.19031935799849300L 05 
0. 412412876 214676 30D 06 
•0.47 843 24 3 277 376 570D 07 
0. 3 3 3826 97661047 380D 08 
■0.14 98 36 5 806 4 213 5 30D 09 
0.44986285138408280D 09 
0.92126 5609095 21140D 09 
0. 1290096 25 69385 440D 10 
•0.12149596007792130D 10 
0.73558445760512930D 09 
■0.25859050142^34330D 09 
0.401165 9999999 99 90D 08 

1 = 8 

•0.22422905954241500D 01 
0.465853157890122800 03 
C.23441235145503930D 05 
0.50244970703015500D 06 
■0.5 74 16 32 37 8045 85 60D 07 
0.39324448510673690D 08 
•0.172817442907787 SOD 09 
0.50733061810099180D 09 
■0.1015 62016 907 3 3 83 OD 10 
0.13911041045129520D 10 
■0 . 128284 39 3 9002 5 3 80D 10 
0.7615932405 37007 40D 09 
■0.262925 2985 7765 640D 09 
0.4011b599999999990D 08 



( N ■ 14 ) 



1 = 9 

-0.29373422408393690D 01 
0 . 608213691132445 50D 03 
-0.30289244915877940D 05 
0.63808107291382380L 06 
-0.7 12299 3035 6 286 800D 07 
0.47456473659044190D 08 
-0.20249980523924240L 09 
0.57740263089002750D 09 
-0.11C46104691522120D 10 
0.1502l6d3094537710D 10 
-0. 135 4 25 85 45 77 5815 OD 10 
0.78789815524134710D 09 
-0.26715875162966600D 09 
0.40116599999999990D 08 

I = 10 

-0.4125 8 26 45 45 6 9 95 80D 01 
0.8494011115 25 44690L 03 
-0.41549543307381450D 05 
0.84979835908124590D 06 
-0.91315334289765970D 07 
0.58371048055375070D 08 
-0. 2 39402036 989 405 OOD 09 
0. 65 87745 181819 2800D 09 
-0. 1243915 1879965760D 10 
0.16177332191612960D 10 
-o. 1425407 1787365220D 10 
0 . 81314437068440010D 09 
-0. 2710929117226626 OD 09 
0.401165 99999999990D 08 
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COEFFICIENTS OF THE POLYNOMIALS 
P n (x)/(x-x i ) (continued) 



1 ( N » 1 4) 

I - 11 

-0.6 3 9576 149941905 OOD 01 
0.13022041497201630D 04 
-0.61513128408473B90D 05 
0.11896 55 5 8795 21140D 07 
-0 . 11981847 811414090D 08 
0.72255528335757830D 08 
-0.28231372004949410D 09 
0. 74676425157665 550D 09 
-0.13655274325948350D 10 
0.17299700043088120D 10 
-0.14917937922113510D 10 
0.8 35 9247986 3 585 380D 09 
-0.27454382727258420D 09 
0.4011o599999999990D 08 

I = 12 

-0.115 7416 2157 5 8615 OD 02 
0 . 2 2966 12 82 34427 870D 04 
-0.99808481529H3370D 05 
0. 17096 37067 3 2516 90D 07 
-0.15 66 473649057 2830D 08 
0.88131683554205400D 08 
-0.32713902258233050D 09 
0.83257505630568550D 09 
-0.147795 3998849806 OD 10 
0 . 182941648212205 20D 10 
-0.15 486 3185 367 37 270D 10 
0 . 8548878292 2267090D 09 
-0.27735015213806370D 09 
0.40116599999999990D 08 



(N g 14) 



1-13 

-0 . 27 9465 765 2 197 5 180D 02 
0. 5 087769 9 30 3 1616 40D 04 
-0.16299086393618340D 06 
0.2 3 6229996 9 3 43 89 30D 07 
-0. 19583843820272490D 08 
0.10327411926099290D 09 
-0. 36 671946 9 2 26 316 20D 09 
0.90416930565976880D 09 
-0. 15 6780306 41774040D 10 
0. 1906 27 27 57 30076 80D 10 
-0.15914037904035320D 10 
0. 86885 374294133100D 09 
-0.27938072542698890D 09 
0.401165 9 999 9999990D 08 

I = 14 

-0.14561307272487430D 03 
0.93592930947540420D 04 
-0. 2 2757146 9476 83 940D 06 
0.29087565 319279450D 07 
-0.22499462811690370D 08 
0. 113704040877 34690D 09 
-0.3925 6 36 5 89841 3480D 09 
0. 94914H8241443470D 09 
-0.16226264769256810D 10 
0.19521312846090240D 10 
-0.1616 47 95 5 2657 9080D 10 
0. 87592672246302460D 09 
-0. 2805 410765 199 8 310D 09 
0.40116 5 9999 9999990D 08 
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COEFFICIENTS OF THE POLYNOMIALS 
P n (x)/(x-x i ) (continued) 



N = 1 5 
1 = 1 

-0.10060401524772970D 01 
0.2404374b422705070D 03 
-0.141243615240087B0D 05 
0.35931285982064320D 06 
-0.496l212977893b520L 07 
0.418485468582409501) 08 
-0.23113042032021540D 09 
0.871552796931459600 09 
-0. 2 297410852 7756 610D 10 
0. 427 22 96 4284119 88OD 10 
-0 . 5 577274 3622 88OOOOD .10 
0 . 4998946 4172 22 b740D 10 
-0.29282905809682000D 10 
0.10091951654130530D 10 
-0 . 15511751999999990D 09 

1 = 2 

-0.10323781804559840D 01 

0 . 24670510910533980D 03 

-0. 14487676 2 8 3525 810D 05 
0.36834483466923620D 06 
-0.50817764644264680D 07 
0 .42619702885 791220D 08 
-0.236178969898214200 09 
0.889158622092550000 09 
-0.2 3 39 3827 3 946 906 90D 10 
0.4 340817 147 46 56 100D 10 

-0.56525518329132970D 10 

0.50521073445626440D 10 

-0.295007299088150300 10 

0.10131288779044060D 10 
-0.15511751999999990D 09 



(N- 15) 



1*3. 

-0.10821295604528690D 01 
0. 25854022227177740D 03 
-0.15173043661639020D 05 
0.38535406088777140D 06 
-0. 5 30826 5843 3782880D 07 
0.44638133558104910D 08 
-0. 24559287741861270D 09 
0 . 9 218209 37 36 457 290D 09 
-0.24167749132397410D 10 
0.44662567959738930D 10 
-0 . 57 89 2115 87 1097 700D 10 
0. 5 1476 895013676 600D 10 
-0.29888074454704940D 10 
0.10200367891668490D 10 
-0.15511751999999990D 09 

1 = 4 

-0.11598116781339040D 01 
0.27700966416469860D 03 
-0. lo 24083 315 9132670D 05 
0. 41177 8608377 9 3960L 06 
-0.5 65 8676 9004 905 5 30D 07 
0.474361105 847 00210D 08 
-0.25997785262680260D 09 
0. 97 13108049 66 756 50D 09 
-0.25328663064844250D 10 
0.46522268248212260D 10 
-0. 5 9891050 3095 5 1950D 10 
0.52853821741328200D 10 
-0.30436640474547060D 10 
0.10296376990336900D 10 
-0.15511751999999990D 09 
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COEFFICIENTS OF THE POLYNOMIALS 
P n (x')/(x-x^) (continued) 



(N = 15) 



1=5 

-0.12730970773184230D 01 
0.30392251000383360D 03 
-0.17792902679660220D 05 
0.45002337213037700D 06 
-0.61627019399481760D 07 
0.5 14277 85 6 277697 50D 08 
-0. 2302895407 200716GD 09 
0.10403253648001450D 10 
-0.26924031654720320D 10 
0.49035 297017 800650D 10 
-0.62541611670470500D 10 
0.5 46 419256 5 3287 200D 10 
-0.31133126064361470D 10 
0.10415 387 46 2979740D 10 
-0.155H751999999990D 09 

1 = 6 

~0. 14345645 27 35 85 880D 01 
0.34223749949777700D 03 
-0.19994618991286290D 05 
0.5 039415 29 36 85 0410D 06 
-0.68669710511210130D 07 
0.56940074737338770D 08 
-0.30793105169337000D 09 
0.1132616756368L360D 10 
-0.29014838853350630D 10 
0.52255067461657440D 10 
-0.65854823344066340D 10 
0 .5681915 922285 2080D 10 
-0. 31958357783530470D 10 
0.10552527502683340D 10 
-0.15511751999999990D 09 



(N = 15 ) 

1 = 7 

-0.16o50098603607240D 01 

0 . 39683010652589540 D 03 

-0.23115614641881100D 05 
0. 5 79 6971315 426 105 OD 06 
-0.78439327865484000D 07 
0.64460155784839010D 08 
-0. 34487542 895 894950D 09 
0.12530451403345920D 10 
-0.3167057o123719550D 10 
0.56227 381145203220D 10 
-0.69819660060510230D 10 
0.5 9344986 9885 3 2170D 10 
-0 . 3 2 886 5 6 86406046 80D 10 
0. 10702 18275 5502 340D 10 
-0.155H751999999990D 09 

1 = 8 

-0 . 20000000000000000D 01 
0 . 47 6 00000000000000 D 03 
-0.27607 9999999999 90D 05 
0.687 34400000000000 D 06 
-0.92067919999999990D 07 
0.74703439999999970D 08 
-0.39377575999999980D 09 
0.14073497599999980D 10 
-0.34956416599999970D 10 
0.6096 8317 199999950D 10 
-0.7438509H99999960D 10 
0.62153985599999960D 10 
-0.33885154799999980D 10 
0.10858226399999980D 10 
-0.15511751999999990D 09 
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COEFFICIENTS OF THE POLINOfelALS 
P (x-)/(x-x i ) (continued) 



(N = 15 )- 
1 = 9 

-0.25037371216725420D 01 
0 . 594628209665652 }0D 03 
-0.34264573377579380D 05 
0.843798034265756001) 06 
-0. 1113 397367 825 6 5 90D 08 
0.88693731620030360D 08 
-0.45792748152110470D 09 
0.16011978721570450D 10 
-0.38907 3917 99828800D 10 
0.66432116c39475960D 10 
-0.79440440592891930D 10 
0.65151354569540330D 10 
-0.34915136019074060D 10 
0. 11014 270044 4976 20D 10 
-0.15511751999999990D 09 

I = 10 

-0. 33011544884820390D 01 
0.7 8137945 6 30398640D 03 
-0.44561031797570630D 05 
0. 107 85 497 88 38834 30D 07 
-0 . 13905090623848970D 08 
0.107 79 3988 5 4 364 1800 09 
-0.54072029595592120D 09 
0. 183785 287 417 4 3480D 10 
-0.43486692904945690D 10 
0.72473157080254900D 10 
-0.84799085038268430D 10 
0.68211751770212770D 10 
-0.35932443448646990D 10 
0 . 11163925^97 316650D 10 
-0.15511751999999990D 09 



(N= 15) 



I = 11 

-0.4oo17022773344^50D 01 
0.10970770784429030D 04 
-0.61454861805635260D 05 
0.14443125523029000D 07 
-0.17930899592615290D 08 
0.13345340595477830D 09 
-0.643957 82849 385780D 09 
0.21140484721628930D 10 
-0.48534013477272510D 10 
0.78813357285103280D 10 
-0.90192310519957370D 10 
0.71183137806629320D 10 
-0.36890032245624610D 10 
0. 11301065 337020220D 10 
-0.15511751999999990D 09 

1 = 1 ? 

-0.72573609683551000D 01 
0. 1689097344182373OD 04 
-0.91376725490805910D 05 
0.20313590993403970E 07 
-0.23654503729479240D 08 
0.16622265564157130D 09 
-0.7 6 46984 320035 4 300D 09 
0.24149028864505400D 10 
-0.53723899378805530D 10 
0 . 85032145783714720D 10 

-O. 9527963843868932 OD 10 

0.7 3896 2 3 15 06 95 6810D 10 
-0.37740682800167090D 10 
0.11420075809663060D 10 
-0.15511751999999990D 09 
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COEFFICIENTS OF THE POLYNOMIALS 

4. 

P (x)/(x-x i ) (continued) 



(N*. 15) 



I » 15 

-0.13175801987562450D 02 
0. 29885 9071899985 30D 04 
-0.14877337284699890D 06 
0.29317032602886430D 07 
-0.31082200963825030D 08 
0.203 9128097 105 89 80D 09 
-0.89171865038657040D 09 
0. 27 106 839576 540920D 10 
-0 .5 85 6 46785 390876 20D 10 
0. 905 9 5 45 2 89113 26 40D 10 
-0.99680342913989280D 10 
0 .76 17 907 5 6 926 393 50D 10 
-0.38440235063014350D 10 
0.11516084908331470D 10 
-0.15511751999999990D 09 

I * 14 

-0.318843960568736300 02 
0.66356 403417 380920D 04 
-0.24373579094520960D 06 
0. 4066670076 26 36 210D 07 
-0.390287302497177bOD 08 
0 . 24008254354781750D 09 
-0. 10046 383076 6887 5 OD 10 
0. 29592651630077490D 10 
-0 .62463262239H8420D 10 
0. 94929825630010400D 10 
-0.103019668802427600 11 
0.77874933494715830D 10 
-0.38950918981242010D 10 
0 .115 85 16402095 5 900D 10 
-0.155H751999999990D 09 



(N «15) 



.1 * 15 

-0.166 56 2815 035 8 32 90D 08 
0.12231904255938610D 05 
-0.34113659259376500D 06 
0 . 502077079235 4 36 50D 07 
-0.44966831198449170D 08 
0. 2651148349447 8620D 09 
-0.10787415823338880D 10 
0.3H62333H5351080D 10 
-0.6 485 5027 10946 16 30D 10 
0 . 97 5 305 38298919 200D 10 
-0. 10498862 328606 5 50D 11 
0. 7 88609493 2177 5700D 10 
-0. 39244477505985050D 10 
0 . 116 245 01145 869460D 10 
-0.15511751999999990D 09 
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APPENDIX G 



ELEMENTS OF THE EXPLICIT INVERSE 
MATRIX WITH DIVISION BY WEIGHTS 



From equations (37) and (45) it is evident that the explicit 

inverse matrix with division by weights is the necessary matrix for 

use in obtaining the temperature function: 

N-l 

a(x. ) - u(-log x ) - j> F (k+1) 

<~0 
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ELEMENTS OP THE EXPLICIT INVERSE 
MATRIX WITH DIVISION BY WEIGHTS 



1 = 1 

0 .67620999 2 275 55 400D 00 
/ -0.73524199845511530L 01 
0 . 12000000000000080D 02 

1 = 2 

-0.14999999999999990D 01 
0.14999999999999990D 02 
-0.14999999999999990D 02 

1 = 3 

0.53237900077244490D 01 
-0.16 6475 80015 4 4 8840D 02 
0.11999999999999990D 02 

N = 4 



1 = 1 

-0.65496972587504130D 00 
0.12395555984240360D 02 
-0.45626858219032430D 02 
0. 426 64577 89 3729 35 OD 02 

I_= 2 

0 . 12290561919335860D 01 
-0.22746685660881600D 02 
0 .76 664299881176 180D 02 
-0 . 5 764191995 4 206000D 02 

1 = 3 

-0.24952495419778680D 01 
0.42343845761147220D 02 
-0, 96 26145 99814416 90D 02 
0.5 764191995 4 205 890D 02 

1 = 4 

0.87783059330621940D 01 
-0.49135573227363170D 02 
0. 82 36 6 87 5 46 215 49 30D 02 
-0.42664577893729010D 02 



N = 5 



1 = 1 

0.64457574561626680D 00 
-0. 186 609712 87045 070D 02 
0.115 78146 2040847 30D 03 
-0.23948231451334420D 03 
0.15481333922979940D 03 

1 = 2 

-0.11196 2 22405 308090D 01 
0. 3 2133165 655 799910D 02 
-0 . 19 3 34776 61776 2 280D 03 
0.37563762838230080D 03 
-0 . 217035 5 6145 201720D 03 

i = 3 

0.18749999999999980D 01 
-0.52499999999999990D 02 
0.28874999999999980D 03 
-0.47249999999999990D 03 
0.23625000000000010D 03 

1 = 4 

-0.37321558320700430D 01 
0.95791727360611900D 02 
-0. 36 86 4824974282 2 30D 03 
0.49250461742576680D 03 
-0. 217035 56145 201670D 03 

1 = 5 

0. 13096091215873580D 02 
-0.11370836617381250D 03 
0.32621455387960370D 03 
-0. 379771042405 84490D 03 
0.15481333922979590D 03 
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ELEMENTS OP THE EXPLICIT INVERSE 
MATRIX WITH DIVISION. BY WEIGHTS 



N = 6 
1 = 1 

-J.63670289144860900D 00 
0.26164498964931050D 02 
-0.24117639244238380D 03 
0.823416513^964 3150D 03 
-0.11597231598027260D 04 
0.57023452625729080D 03 

1 = 2 

0.10633062297523300D 01 
-0 . 43 37870246 1460870D 02 
0.39436316819926010D 03 
-0. 13115 64025 46 2 5750D 04 
0. 17703674113121060D 04 
-0.81606492177331800D 03 

1 = 3 

-0.16 2092 92 95 3026 820D 01 
0.65461713775432130D 02 
-0 .575089 1902 2082130D 03 
0.179456 3 95 767 37070D 04 
-0. 22082425 39 2646 3 50D 04 
0.92756392549416520D 03 

1 = 4 

0.26369381625442730D 01 
-0.10382467676719200D 03 
0.83478670215413420D 03 
-0.22372330555568110D 04 
0.24295770882061870D 04 
-0. 92756392549416430D 03 

1 = 5 

-0. 5 21376 395 6 2 34085 OD 01 
0.18819940606882500D 03 
-0.1w787736580489900D 04 
0.23907435979473020D 04 
-0.23099571975544600D 04 
0.81606492177330980D 03 



(N = 6) 



1 = 6 

0.18277283962093760D 02 
-0. 22o 341247 84499930D 03 
O.973O7O453O034343OD 03 
-0.16868691372583870D 04 
0.16914414714836860D 04 
-0.57023452625727700D 03 

N = 7 

1 = 1 

C. 6 350567186 2000200D 00 
-0.349H537906429760D 02 
0.444 279785 35 2146 90D 03 
-0.22113580870536630D 04 
0.5 065 80745 4 875 5 810D 04 
-0.5 364185 369 2 3 8286 OD 04 
0.21240546327832390D 04 

1 = 2 

-0.10301560713220060D 01 
0.56505693626735990D 02 
-0.71390602118058990D 03 
0.35067953361694150D 04 
-0.78710473381123000D 04 
0.80943293742454680D 04 
-0.30785879539965880D 04 

1 = 3 

0. 148785 6424026 55 40D 01 
-0.81203287902126410D 02 
0.10092970779827400D 04 
-O.48131388752HO05OD 04 
0.103 37417 3382135 40D 05 
-0.10039689543210950D 05 
0. 35 8934988926 27010D 04 
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ELEMENTS OF THE EXPLICIT INVERSE 
MATRIX .WITH DIVISION BY WEIGHTS 



(N = ?) 



1-4 

0.21874999999999980D 01 
0 . 11812500000000000D 03 
■0.1417500o000000000D 04 
0. 63524999999999990D 04 
•0.12560624999999990D 05 
0.11261249999999980D 05 
•0.37537500000000020D 04 

•1 = 5 

0. 35204555588710010D 01 
•0. 185 295 214805775 40D 03 
0.20377 37 38846 175 30D 04 
•0 .79266 32 8307875080D 04 
0.13979217961099200D 05 
•0.11496409792365190D 05 
0.35893498892626960D 04 

1 = 6 

•0.69410653211762830D 01 
0.33499054543956470D 03 
■0 . 2 65 5 3296088702 320D 04 
0.86058593537967680D 04 
•0. 1357821977686 3660D 05 
0.10377198349745970D 05 
•0. 307 85 87 95 39985 6 80D 04 

1 = 7 

0.24321935531217310D 02 
■0 .406 2045 416 472 346 OD 03 
0.24240160528103890D 04 
■0.68911106957305150D 04 
0.10105700100432590D 05 
•0.73801424274610320D 04 
0.21240546327832080D 04 



N -.8 . 

I - 1 

-0.63263629050388220D 00 
0.44904361134180810D 02 
-0.75 13077246 126 4850D 03 
0.50790321475309150D 04 
-0.16738927934245450D 05 
0. 2 85 1734961913 2120D 05 
-0. 240995 5625 1708940D 05 
0 .79803 68487 595 4930D 04 

1 = 2 

0. 10088738950655580D 01 
-0. 715 15 86 960827 35 20D 02 
0. 11915 7159704 86O8OD 04 
-0.7 995 5 6 9867 3191 210D 04 
0.26057030694553270D 05 
-0 . 43705585221165210D 05 
0.36178290078771380D 05 
-0.11664144752831990D 05 

1*3 

-0.14081395444254550D 01 

0.995399513236I8810D 02 

-0.16437571873o41720D 04 
0 . 1085 6 21461398125 OD 05 
-0. 34559345638578920D 05 
0.56179519216495980D 05 
-0. 4 47 4966106 3421930D 05 
0. 1382342576 9235650D 05 

I a 4 

0 . 19506 2715 411025 6 OD 01 
-0.13714860269817760D 03 
0 . 22260096024338430D 04 
-0.14261847226419680D 05 
0 .4348676 385 994 3770D 05 
-0.O7093136369321480D 05 
0.50629390880600230D 05 
-0. 14854809783306000D 05 
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ELEMENTS OP THE EXPLICIT INVERSE 
MATRIX WITH DIVISION B£ WEIGHTS 



(N = 8) 

1 = 5 

-0. 2 827011b 1340087 10D 01 
0 .196 6206 8346 25 1010D 03 
-0. 30804548507999980D 04 
0.185 76 6 802 835 6496 OD 05 
-0.52456397193370430D 05 
0.7 5 2677 96 5 35 146 170D 05 
-0.53354277602541830D 05 
0.14854809783305980D 05 

1 = 6 

0.452752212711277102 01 
-0 . 306 896 9508214 87 60D 03 
0.44H0301906861850D 04 
-0.2 3 24070487 94 35 370D 05 
0.58913236415b20400D 05 
-0.7797 3493 98991 3020D 05 
G.52014319321227460D 05 
-0.13823425769235560D 05 

1 = 7 

-0 . 891446 665 6 27 893 40D 01 
0. 5 5 415 8402 277 104 60D 03 
-0.57814946771425290D 04 
0 . 25 5025 6 40744 5 03 90D 05 
-0. 58041610578821550L 05 
0 .715 828045 5 6006030D 05 
-0.45470723191052030D 05 
0. 116 641447 5 2631 820D 05 

1 = 8 

0.31230068535155660D 02 
-0.675663616O9512550D 03 
0 . 5 32011176 76 994600D 04 
-0.20618588633532940D 05 
0.43667373451622040D 05 
-0.51507750348382530D 05 
0. 3176 302 316 145 87 40D 05 
-0.79603684875953020D 04 



N = 9 
1 = 1 

0.6 309469015 9499 3 50D 00 
-0.5 6 144067 147 71 89 30D 02 
0.1192^225317072080D 04 
-0.10446389121829060D C5 

0 . 402245892342210300 05 

-0.1121O523386107140D 06 
0.15126259454114080D 06 
-0.106139538456527600 06 
0.30l8o269946095430D 05 

1 = 2 

-0.994349226306987500 CO 
0. 8o4O027 96 2 o242 2700 02 
-0.187250776427494800 04 
0.16335839433292440D 05 
-0. 717861912321406 30D 05 
0.17 26 2 94405 3493 840D 06 
-O.22y99798630OO793OD 06 
O.15897455931633O70O 06 
-0 . 443817G005984H40D 05 

1 = 3 

0.135612170050724700 01 

-0.120369850651937700 03 
0 . 25359056032479280D 04 
-0.219175186793261400 05 
0.950031686422515100 05 
-0 . 22431467127417700D 06 
0 . 29207121661956960D 06 
-0. 196 4419 P 8 j 2 8 307000 06 
0.531885302580274300 05 

1 = 4 

-0.180760131953163100 01 
0.159954125710033600 03 
-0.333747430577061100 04 
0.28363935208562880D 05 
-0.120009178736693300 06 
0.27472299293684350D 06 
-0.345041748289709900 06 
0.22333H7187123050D 06 
-0. 56191387545852760D 05 
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ELEMENTS OF THE EXPLICIT INVERSE 



MATRIX -WITH DIVISION. BY WEIGHTS 



( N_- 9 ) 

I =. 5 

0.24609 374 999999980D 01 
-0.21656250000000000D 03 
0.44395312499999990D 04 
-0.36599062500000000D 05 
0.14850773437499990D 06 
-0.32376093749999990D 06 
0.387105 46874999970D 06 
-C.23930156249999990D 06 
0.59825 390625000010D 05 

I *= fi 

-0.35423371992380150D 01 
0.30832612996597780D 03 
-0.61012780666898300D 04 
0.47404503239454910D 05 
-0 .1788265 5 111494290D 06 
0.36429059007332490D 06 
-0.41108239647497120D 06 
0.24219992849559070D 06 
-0.58191387545852670D 05 

1 = 7 

0 .5 65 89 9214085 4784 OD 01 
-0.48003575 926795 OOOD 03 
0. 87216161808651660D 04 
-C. 5 946 192 1498086 270D 05 
0.20222u63113608010D 06 
-0.38138919379832560D 06 
0.40625635414618650D 06 
-0.2290662835359H80D 06 
0.53188530258027340D 05 

1*8 

-0.11134161372277090D 02 
0.8662o 6 307 803 2 3180D 03 
-0.11479412061649910D 05 
0.657 39913426 9945 70D 05 
-0 . 201218241175 93400D 06 
0. 35 4267946 97368490D 06 
-0. 35986367476130570D 06 
0.19607904116239290D 06 
-0. 4438170005 9839920D 05 



(N = 9 ) 



1 = 9 

0.390016 942 8945 2010D 02 
-0.10602788886162110D 04 
0.10622421003920810D 05 
-0.5 3 50879107 4806 880D 05 
0.15253238835015280D 06 
-0. 25 7003142 8247 9870D 06 
0.25355738385848900D 06 
-0.13535662111862080D 06 
0. 3018 826 9946 892 2 30D 05 

N - 10 



1 = 1 

-0.6 297 207 88895 116 80D 00 
0.68631241594361980D 02 
-0.18007322514163060D 04 
0.19787480981637290D 05 
-0.11232455134290630D 06 
0.36273559299210630D 06 
-0. 6 9145 81879495 07 20D 06 
0. 76 901847 43085 2180D 06 
-0.46080521705457350D 06 
0. 114 826775 61427 35 OD 06 

i = 2: 

0.98397308015031300D 00 

-0.107l8l875o6752850D 03 
0.2 80746 36043 399 86OD 04 
-0. 30759373574475980D 05 
0.17385617851724700D 06 
-0.5 5 8192 9246726 4 940D 06 
0. 105 6 1498767 43 17 30D 07 
-0.11637950926080640D 07 
0. 6895 5 60021286 2 64OD 06 
-0.16952953248708180D 06 



^ LEI. ENTS u? I HE EXPLICIT 1 IN ; SHSE 
Iw-.TKIX WITH DIVISION PY WEI HITS 



( N = 1 0) 

I = 3 

-0 . 13 2009-^2 620004000D 01 
0 . 14 36 3 80600905 6600D 03 
-0 . 37496 16 269 3 ob 5 9 10D 04 
0.40840168676119880D 05 
-0.228560256d3o27430D 06 
0 .726 4 3928184 4946 50D 06 
-0.135486007841809400 07 
0.146728284050193300 07 
-0 . 85 20275 15 E 3799080D 06 
0.204799772 «u030410D 06 

1 = 4 

0 . 171 3486578 5 o5H50D 01 
-0. 186094715 33044610D 03 
0 . 4 82 9402 1180 3102 10D 04 
-0.5^0684504051792700 05 
0.2o75415069b^2o020D 06 
-0.89548791099262010D 06 
0 . 16 3 207 2 1886 5 4026 OD 07 
-0.17216549991630^700 07 
0.97183826091366580D 06 
-0 . 22688995 36 34 1056 OD 06 

1=5 

-0 . 22391285740 3962 10D 01 
0.242406191036718700 03 
-0.822822282254057000 04 
0 . 06OO4 5 67 96 10048100 05 
-0. 3 5 378415 85 76985 800 06 
0.107511267183415100 07 
-0 . 189 39045 89 842 8100D 07 
0.192861439211779600 07 
-0 . 105 169 2 1388897 3 90D 07 
0.237640313597974800 06 



(N = 10) 



1 = 6 

0.30^244131296665000 01 
-O025366323329349IOD 03 
0.821209538592971400 04 
-0.844331636655060400 05 
^.436943888137282300 06 
-0.12605^ 710277020500 07 
0 . 2 12080260829802 oOD 07 
-0.2^701^57052o9370D 0; 

0 . 10870706 034920^600 07 
-0.237o4031359797450D 06 

1 = 7 

-w .4334775502 c589230D 01 
O. 40I3 244 2 191695 8900 03 
-0 . 112 45 195 102 879 3 30D 05 
0.109076^24196158500 06 
-0. 5 26 19608282 311960D 06 
0.142300130821287400 07 
-0.22o67976051638040D 07 
0 . 21149 87 24267 9 9040D 07 
-0.10701713^179308300 07 
0.226889953634105608 06 

1 = 8 

0. o915289 3 65 90 37 2900 01 
-0.71754087313121890D 03 
O.I 0O6 20 38367213 5 80D 05 
-0.137129876124356000 06 
0.596179707423714900 06 
-0.150144698894784400 07 
0.22625 302986972490D 07 
-0.202385454156894600 07 
0.99117^440504742600 06 
-0.204799772860302800 06 
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ELEMENTS OP THE EXPLICIT INVERSE 
MATRIX WITH DIVISION BY WEIGHTS 



I N = 1 0) 

1 = 9 

-0. 1 360025 147 207 96 10D 02 
0.12944478368158440D 04 
-0.21206736297286420D 05 
0.15243922236262340D 06 
-0. 5 9948947 9 34917 960D 06 
0.1406575 58315 22120D 07 
-0. 2023 3 28440826 57 90D 07 
0.17504102451139530D 07 
-0. 83620979025 508760D 06 
0.16952953248707710D 06 

I = 10 

0.47636818940463100D 02 
-0.15888058270619700D 04 
0.19703318379258890D 05 
-0.12468504388204750D 06 
0. 45 6 9 35 7287 5 15 5 4 30D 06 
-0 . 1026 45 5 998116 2 400D 07 
0.14345742062o08630D 07 
-0.1216 3406 5 9 985 6000D 07 
0.57263576347380240D 06 
-0.114826775 61425 5 SOD 06 

N = ,1 1 

1 = 1 

0.62880260230913370D 00 
-0.8236622^639663810D 02 
0.26142904642913460D 04 
-0. 3 5 1228223 308 3765 OD 05 
0. 2477 3476 66 349885 OD 06 
-0.10184725075194570D 07 
0. 25656309060809040D 07 
-0.40097 677 38 36 5 6480D 07 
0. 37879189896732840D 07 
-0.19791458507143100D 07 
O.4387488389408913OD 06 



( N - 1 1) 

1 = 2 

- 0 . 97629228872923490 D 00 

0.12783586044824230D 03 
-0.405 2807 305 67 44170D 04 
0.54340754123511910D 05 
-0.38217790816757760D 06 
0. 15 6 5 12287115 336 80D 07 
-0. 39233653943289440D 07 
0.60947715541134540D 07 
-0.57158374375782260D 07 
0. 29608725820561420D 07 
-0.649O1738642927510D 06 

I = 3 

0.12940094492703080D 01 
-0.16931341404173450D 03 
0. 5 35 55796 82999507 OD 04 
-0.71527329877500880D 05 
0.50020326653883290D 06 
-0.2033H29052140890D 07 
0. 504 85 60312 39 3995 OD 07 
-0 .77 5 36 105 7912 1175 OD 07 
0.71747065316893470D 07 
-0 . 36600707958978520D 07 
0.78967223585514890D 06 

1 = 4 

-0.16477565202582030D 01 
0.21533446996204010D 03 
-0.67853720573359590D 04 
0.90030823198171330D 05 
-0.62369982976702370D 06 
0.25040530553528440D 07 
-0. 61246 301487788360D 07 
0.9241088 289 3189 820D 07 
-0. 8 3826 5 94 3261385 20D 07 
0. 4185 267 3 378896 5 5 OD 07 
-0. 88288361424753310D 06 
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iiLUMTT. OP THE EXPLICIT II; 7 ERSE 
MATRIX. »VI TH DIVISION BY A'ElCtiTS 



( N = 11) 

1 = 5 

0 . 209253200 ^ 9004 940D 01 
-0.2729177129bB;)9 300D 03 
0.65470159711^1c340D 04 
-0 . 112<il27bl08008450D 06 
0.76 5 80444 244 5 34 ^30D 06 
-0 . 3 -/lb 3 37 84 5 4 905 HOD 07 
0.72126476677351980D 07 
-0.10613016080694200D 08 
0. 9 376 603 96 35 30b 5 40D j 1 
-0.4 5 587 73066 836 3 300D 07 
0 . 9 3701110b 44 bb6970D 06 

1 = 6 

-0. 27070312 Ov>O000w30D 01 
0. 3 519140u 24 999 9 990D 03 
-0.10909335937500C10D 05 
0 . 1407656 25 OOOOOOOOD 06 
-0 . 9 3785 09765 6250140D 06 
0.35871304218750020D 07 
-0.830376421875001401: 07 
0. 11821497187 5 00010D 08 
-0.1011647 3554607 500D 08 
0 . 47740661716750060D 07 
-0.95481323437500150D 06 

1=7 

0 . 3636o46o69o 400980D 01 
-0.47010605445151210D 03 
0.14314916797832990D 05 
-0.179234601267..7820D 06 
0 . 1147 471006 2 920090D 07 
-0 . 4 197 42 49 47 24104 80D 07 
0 . 9301840841262 8020D 07 
-0.12725317995043420D 08 
0.109l314bl52l03770D 06 
-0.46113379976303560D 07 
O. 937 OIIIO 64466602 OD 06 



( N = 1 1) 

1 = 8 

-U. 52O49914O3482597O0 01 
0. 6b 541^17 34 347627 JD 03 
-0 . 1956207 36195o7 3oOD 05 
0.23125634510221690D 06 
-0 . 1302631 35 937 56070D 07 
0. 47 5 27 5820427 66 3 10D 07 
-0.999457864598442300 07 
0. 1 3096 5 9671726 8040D 06 
-0.10445016032745860D 08 
o.4o435bo804585649JD 07 
- 0 . 0020636142475^9100 06 

1-9 

0.62vbo42lll262ol90D 01 
-u.1o33oo 59056676870D 04 
0.2793150996854O430D 05 
-0.29127988457733320D 06 
0.157637962104320000 o7 
-0.5 04447 57 345765650D 07 
0.10C50291620006940D 06 
-0.12642161324688600D 05 
0.97693199820902740D 07 
-O. 42366 5 156 265 35 970D 07 
0.759672^3585514^000 06 

I = 10 

-0 . 16 3127 95060626 460D 02 
0 .166 4406 876 5 925 46 OD 04 
-0 . 3b96524o44 5 1012 30D 05 
0.3^513164^79571ol0D 06 
-0.159036 o5149134350D 07 
0.47557994039214340D 07 
-0 . 905 17 570 2 5 15 64 250D 07 
0 .1101o 601364003760D 08 
-0 . 53097o65 85 3 900020D 07 
0. 3 5 37 301282 2 36 47 80D 07 
-0 . 64 ^ 81736642925 160D 06 
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ELEMENTS OF THE EXPLICIT INVERSE 
MATRIX WITH DIVISION BY WEIGHTS 



( N = 1 1) 

I * 11 

0.5 7 135 446 067 44 3050D 02 
-0 . 2^9 319495 037 1555 OD 04 
0.34449271098960500D 05 
-0. 26691108244 36 8290D 06 
0.12171732421673420D 07 
-0.34859840680922630D 07 
0.6 4479931659575 390D 07 
-0. 76 94 194226 2113 810D 07 
0.57193040855838490D 07 
-0.24083425386942900D 07 
0 .43674863 89408 3490D 06 

N = 1 2 

I = 1 

-0.62809712046008180D 00 
0.97349208725481830D 02 
-0. 367 409620818045 10D 04 
0 .5916423620406 3700D 05 
-0 .5 061 37907041945 80D 06 
0.25 67 390897 5 000990D 07 
-0. 818255 3 6131982310D 07 
0.16 806 96 147 815 9360D 08 
-0.22201731817903030D 08 
0. 18207 31946182 3990D 08 
-0, 84295 81077 5 42 3920D 07 
Q, 16 828 13 214672496 OD 07 

1 = 2 

0.97044216519729790D 00 
-0.15 036 96 686187 1660D 03 
0.56705340420321710D 04 
-0.91185184277295630D 05 
0. 778494488935 5 3960D 06 
-0.39383401739560010D 07 
0.125 0946 925 5 8186 80D 08 
-0.25586333143220000D 08 
0.3363522C804863140D 08 
-0 . 27424315890346060D 08 
0.12611867398301090D 08 
-0.2498417 96 27 5 99490D 07 



( N = 1 2) 

'1 = 3 ' 

-0. 127 44705 976447 86OD 01 
0.197 37725 41105 3 610D 03 
-0.74314329933348770D 04 
0.11917694554194690D 06 
-0.10134999401397680D 07 
0 . 5 1007 87 2028815 800D 07 
-0.16097249175721700D 08 
0.32670693022309940D 08 
-0.42551737919709900D 08 
0.34329427598729830D 08 
-0.15 60021996 016 8470D 08 
0.3049867 3596812220D 07 

1 = 4 

0.15997929592334350D 01 
-0.24755198318340960D 03 
0 . 929644422619035 90D 04 
-0.14842587620106170D 06 
0.12542388006541950D 07 
-0.62600943174096300D 07 
0 . 1955 38286 8024145 OD 08 
-0.39205841039798720D 08 
0. 5 035 64795 48145 220D 08 
-0.40001461744518260D 08 
0.1787565925166146OD 08 
-0. 3433439948167 9000D 07 

1 = 5 

-0.19896588b93697860D 01 
0.30747756727537110D 03 
-0.11500307185105370D 05 
0. 18234946701209440D 06 
-0.15258580461689860D 07 
0.75200497695791070D 07 
-0.23133310008012750D 08 
0.45 57711045 9247 880D 08 
-0.57424066378613350D 08 
0. 4469683942395 3880D 08 
-0.19561620276687970D 08 
0. 3 6797 047140185 5 70D 07 
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Ei^r ol' 1 me explk i 2 iaveioe 

i. a. ’. ax .viTii division ;-y 22 



(N = 12 ) 

1-6 

0.2498112036/ 529 350D 01 
-0.3652653109319b790D 03 
u . 14 3188868b 906 215 uD ub 
-0. 2 ^ 46 104 9^373 2141J0 06 
0.185132441389726801) 07 
-0 . 895242ob939160230D 07 
0.2693^3316142588000 08 
-0.518124627894654100 06 
0.6 3678212531130^4'^.D 08 
-0.483571868594873380 08 
0.206655146350393000 08 
-0.3800635988l2573o0D 07 

1 = 7 

-0.321338l8o90451.i2CD 01 
0. 49 3940742022 827 90 D 03 
-0.181702629695584100 05 
0.280116415721284100 06 
-0 .22544986li6o09910D 07 
0.10593935323414500D 08 
-0.30898406193109330D 08 
0.575936093587351200 08 
-0.686932486602302100 06 
0.507370198560090500 08 
-0. 2114 14 812 34 343^ 90D 08 
0.3O006359881257360D 07 

1 = 8 

0.43050516901156380D 01 
-0 .65796 8ll3u4 9961100 03 
0.237745178688517400 05 
-0 . 35 5719911o2c5 88500 06 
0. 275302519211. 5250D 07 
-0 . 1238862773^7401200 08 
0.346504359897684900 08 
-0 . 6 217 892111c 2444 10D 08 
0.717258>379907 30800 06 
-0. 51464395926094200D 08 
0.209151315775160100 08 
-0.367970471401853700 07 



(N = 12) 

1 = 9 

-0.6193)5/31-44920200 01 
0 . 9301 J2 4427 7o069 30D 03 
-0.3244946575x8243200 05 
0.45 .0097275 3631C60D 06 
-0.33205904613610080D o7 
0.14064233^19780170D 08 
-0.373 ,9230994475820D 08 
0.b4 3627001533227900 C 8 
- 0 . 7 17 o960127 5 4 55 3 30 D 08 
0 . 50O6406637713758OO 08 
- 0 . 13892 ] 901781652600 08 
0.34334 399481678820D 07 

I = 10 

0.9803194840362150oL 0, 
-0.14440o92792927150D 04 
0.4632u0015746l2040D 05 
-0.57863539694709920D 06 
0.380221433697761400 07 
-0.149956697635312300 08 
0. 37 61336705 92 19770D 06 
-0.62546016093489820 L 06 
0.876313266086760300 06 
-0.460699327795117900 08 

0 . 17940320996324770.) 06 
-0.30498673596811690D 07 

I = 11 

-0 . lb 27 182 7 27 6777 '-500 02 
0 . 2 604 417o 5 8 169 901^0 04 
-0.614215396868876100 05 
0.6479296^ r 05o3506600 06 
-0.384694746725385600 07 
0.142084721814040300 08 
-0 . 342 4 2 16 015 lb 5 5 2400 06 
0.54635778o71493580D 08 
-0.578885531401082100 08 
0.38716629659131920D Oc 
-C . 148707 3019205 7 450D 08 
0.24964179627597960D 07 
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ELEMENTS OF THE EXPLICIT INVERSE 
MATRIX WITH DIVISION BY WEIGHTS 

(N = 13) 



( N = 12 ) 

I = 12 

0.67 497 577 8702485 50D 02 
■0.32085918266690160D 04 
0. 57375030486011780D 05 
■0.53340516H7338810D 06 
0 . 295 35'- 36126 9913 lGL 07 
■0.10448804284924580D 08 
0. 24480199520638750D 03 
0.38435^39097431290D 08 
0.39997175270061230D 08 
•0.26466235493383090D 08 
0.10081 36 4283 85 3400D 08 
■0 . 16 82 613214 67 22 190D 07 

N =_ 13 

1 = 1 

0.62754331553197030D 00 
■0.11358033610843180D 08 
0.50^50939417970050D 04 
•0.95442184513324400D 05 
0.97168744238521400D 06 
■0 .5 940496 8247014960D 07 
0.23229631009714890D 08 
0.60063715334341270L 08 
0.10380582544494940D 09 
■0.1185 5 57 87 99997 500D 09 
0.65 633117 6 4 25 66 710D 08 
0 . 35 6 6486 2 301875 86 OD 08 
0.64752096878575660D 07 

1 = 2 

0. 96 5 8807206915 9760D 00 
0.17478290345705200D 03 
0. 77 28269474 ^ 5 6 325 OD 04 
0 . 1466 35 0916 9035 800D 06 
0.14907030198521330D 07 
0.90960216544733470D 07 
0.35483079564821880D 08 
0.91477792834819540D 08 
0.15754736564692130D 09 
0 .17920322086746790D 09 
0.12913513896815610D 09 
0. 5 337199515 8417440D 08 
0. 96 3184643 20360660D 07 



I » 3 

0 . 12594266332923030 D 01 

-0.2 27817 86 91145 6130D 03 
0 . 10061611146213 320D 05 
-0.19054009023605280D 06 
0.19 31643916 25 6 4420D 07 
-0.11743348528528470D 08 
0 . 45600366 924790600D 08 
-0.11691192072352450D 09 
0 . 20004 5 4 67 17 896 5 5 OD 09 
-0. 22 5 6444212 917 45 20D 09 
0.16137022920607280D 09 
-0 . 6606 6 567 924 32045 OD 08 
0 .117992675145 86 200D 06 

1 = 4 

-0.15 6 35 9 30 3788 26 840D 01 
0.26266984060442140D 03 
-0 . 12461600812048860D 05 
0 . 23524972654785430D 06 
-0.23742867586914440D 07 
0.14 3504 35 7 368 8825 OD 08 
-0.55323068110608670D 08 
0.140625255C6500070D 09 
-0.23o24303520694750D 09 
0.26597547632329040D 09 
-0.18771579270607220D 09 
0.75836132834923010D 08 
-0.13354193591661380D 08 

1 = 5 

0. 1914 224o6 475 90760D 01 
-0.34574563103880170D 03 
0.15 20011 2938246 720D 05 
-0 . 26559473259563270D 06 
0.26631024340339590D 07 
-0.17154969863530360D 08 
0.6 5 436 80261295 3 100D 08 
-0. 164 28685 6 2595 0120D 09 
0.27446106194385310D 09 
-0.30181562617390290D 09 
0.209 d0921505548620D 09 
-0.83281068929590510D 08 
0.1441908265 902 9250D 08 
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ii’LjL ,, 0/ i’Kti jSXPi-ICir xivV/. 

.<.;.i\ttX ,VITH DIVISION t:Y -VSIOOPS' 



( N = 13) 

1 = 6 

-0 .235110206 376119400 01 
0.424079069301046000 03 
-0.165662232690400000 05 
0.34637477026u33250D 06 
-0.34378691560332^900 07 
0.203376910794690200 08 
-0 .76 4067 1482 2997720D 08 
0.188561468469960600 09 
-0.30924414061361730D 09 
O.33353320311064930D 09 
-0.22715895046998^900 09 
0. 005 31198296 3 111300 08 
-0.156441199260217500 06 

1 = 7 

0.^32bl71o75000020D 01 
-0 . 52767109375000080D 03 
O.22962392570125O2OD 05 
-0 .4^3763l8359375070D 06 
0.414290830070l25o0D 07 
-0.240522000466749800 08 
0.664338924216701600 08 
-0.213241622343750300 09 
0.341544650712691200 09 
-0.3599115441796681^0 09 
0.239737689597656600 09 
-0.915029349609376600 08 
0 . 15 2 5048916015 62 oOD 06 

1 = 8 

-0.3759293546021u4000 01 
0 . 67442121115 956970D 03 
-0.290358274456996100 05 
0.5266^500414 O-910D 06 
-0 . 5 O'. 65 15939c 6 17 1000 07 
0.2o38493052929ol20D 08 
-0.101256110701616800 09 
0. 36 91794 o05 65 30400 09 
-0.366789729244647200 09 
0 . 37 85 4677 945665 700D 09 
-0.246227664459996900 09 
0.9199824O639949890O 08 
-0.150441199260217200 08 



( N = 1 3) 
l = 9 

0.5 t / 5 d0o 326371340D 01 
-O.89u7o5 387598239000 03 
0.3792 300I62396146 OO 05 
-0.667688055269478100 06 
0.613409249322276700 07 
-0.331942497720542500 08 
0.1136980 37741227400 09 
-0.256465o22cs7570360D 09 
0 . 3 600 24 6 467 127 40100 09 
-0.3Ob01>91623^93430D ^9 
0.^45176912325921700 09 
-0.897479229787600400 08 
0. 14419' ''826590291600 08 

I = 10 

-0 .7100095 765 o2c5 0000 01 
0.1266b^5985796 J150D 04 
-0.517219499684607100 05 
O.Oo0732465039ol200D Ob 
-0. 7405 1109220O 630700 07 
0.377652683642^23900 06 

-o . 1^3109234714406300 09 

0 . 2b 6b 94 49291O41760D 09 
-0. 3690 3833027967030D 09 
0.378117735026167300 09 
-0.234695108564767900 09 
0.844141902674127700 08 
-0.1336 1935918612300 06 

I = 11 

0.1143502^541 3891900 02 
-0. 196 5 914425 9 8146 60 1, 04 
0.73637337089061/ 300 05 
-0. 1087 167402 <i/405 900 07 
0.650063755954654700 07 
-0.404117275349861200 06 
0.12,066330081462900 09 
-0.2605^.6011460648700 09 
0 . 3 6 67 5 o 4020 3 3629 oOL 09 
-0.350034386140314900 09 
0.213369418001231200 09 
-0.755246222507117900 06 
0.11799267514565840D 08 
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SLANTS OF THE EXPLICIT INVERSE 
MATRIX WJTH DIVISION BY WEIGHTS 

N - 14 



(N = 13 ) 

I at 12 

0.22477370421262190D 02 
0. 35 45 3247719 12 66 10D 04 
0.9803976 32 996 2 95 60D 05 
0.12 20415 96 68198080D o 7 
0.86285698775474890D 07 
0.38431084088930520D 08 
0. 11 37 37738292 822 80D 09 
U.22946562918523630D 09 
0.31718441412569740D 09 
0.2 95 6 946 5 014905 7 60D 09 
0.17774505 67 399 3&10D 09 
0. 6221016 2026013170D 08 
0.96 318464 3 20357 3 50D 07 

1-13 

0.78723215761735440D 02 
■0.43733370845259480D 04 
0.9 1749 24266 5 078710D 05 
•0.1006 9 24 3767794 210D 07 
0.66410657217231400D 07 
■0.28341539376072510D 08 
0.81542617182726970D 08 
•C.16131015049626300D 09 
0. 2198205 430405 8640D 09 
•0. 20^7 5 409 315 112 2 90D 09 
0.12088347172049520D 09 
•0.42037653952402330D 08 
0.6475 2096o7 65 5 6 300D 07 



1 = 1 

-0. 6 271005 8676 8491 30D 00 
0.13105968903823980D C3 
-0. 67159735263579190D 04 
0 . 14645 75 83114176 80D 06 
-0.17713639267792060D 07 
0 . 12B14838745636220D 08 
-0.6 0088842 69786 9 890D 08 
0.16975509790027530D 09 
-0.41111999022164400D 09 
0.61198659250115030D 09 
-0.61492202573029410D 09 
0.39829832256181480D 09 
-0.15007902306178370D 09 
0 . 249646 1269 967 85 30D 08 

I = 2 

0. 96225 35 6315 207790D 00 
-0.20107528446341820D 03 
0.10299^71598740270D 05 
-0.22749551893141230D 06 
0 . 2711502407462771CD 07 
-0. 19 5 884 19447 5 3 89 5 OD 08 
0.91687380045722230D 08 
-0 . 2 889194 35 2 5 2649 80D 09 
0.6 24 381 I 925 7 867 880D 09 
-0.92670839285596560D 09 
0 . 928013821860426 30D 09 
-0.5 98797396646 43370D 09 
0.22465816840735830D 09 
-0.37221050666135830D 08 
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ELEMENTS OP THE E.P..TCXT I N7ERJE 
R.ATRIX WITH DIVISION 13/ iVEIHHPS 



(N=l4) 



1 = 3 

-0.12475865212784120D 01 
0.260627607177d 7190D 03 
-0.13338370066251030D 05 
0.294202b443612^140D 06 
-0.34994068475626250D 07 
0.25212540726810810D 08 
-0.1176l7540bl302880D 09 
0.36913741442143300D 09 
-0.79397229603140680D 09 
0.11720111890433050L 10 
-0.11664350161564310D 10 
0.74745627373425050D 09 
-0 . 27b 399005 ^96 29120D 09 
0 .457247 3635 6779090D 08 

1 = 4 

0 . 15 35 5 2882846045 30D 01 
-0.32064094O43427940D 03 
0 . 16 3 879094o 21 15 83 OD 05 
-0.36064900873599210D 06 
0.42759286645134930D 07 
-0.30677579805035340D 08 
0. 142 36 676 9 906 5125 OD 09 
-0.44403634576476030D 09 
0. 9481913346704365 OD 09 
-0.13882270207062360D 10 
0 .13690686886 524 820D 10 
-0. 86658664432220080D 09 
0 . 319 9 3822789747 5 40D 09 
-0.5196878727 3840710D 08 



(N = 14 ) 

1 = 5 

-0.18570160912215590D 01 
0. 38752227419995950D 03 
-0. 19767 119057 866 770D 05 
0.4 3 35 576 97 1903 9980D 06 
-0.5 115 8922 177 818850D 07 
0.36477406002075480D 08 
-0.16800267212292770D 09 
0. 5 193355 9 36 38 380 30D 09 
-0. 1097756 5 1237 369 30D 10 
0.15891593037308090D 10 
-0. 15 481701995 747 85 OD 10 
0.96954303861109500D 09 
-0.35232510694594320D 09 
0 . 5 644086 881409 36 20D 08 

1 = 6 

0 . <i' 242976 7477604 370D 01 
-0 . 46762438068047850D 03 
O.23784307536616040D 05 
-0.51912052591397640D 06 
0.6u83296839682C860D 07 
-0.42991511298815780D 08 
0. 195 8916 64 3129 2110D 09 
-0 .5 98106097 57789 330B 09 
0.124703057774ol890D 10 
-0.17788305398299360D 10 
0 . 17064 3926 9870c ObOD 10 
-0 .10519615804934060D 10 
0. 3762879795666 96 40D 09 
-0 . 59 347 2618816799 30D 08 
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ELEMENTS OF THE EXPLICIT INVERSE 
MATRIX WITH DIVISION BY WEIGHTS 



(N « 14 ) 

1=7 

-0. 27347513435 8l8'j4^D 01 
0.5 6936165 211156270D 03 
-O.288358O6993646890L 05 
0.6 24858039943 37690D 06 
-0.72488607759361210D 07 
0 .5057903918995 5630D 08 
-0.22702150561159070D 09 
0.6 815 995 226 4 16 5 800D 09 
-0 .1395 834747066 85 30D 10 
0.19546601596146270D 10 
-0. 184081 855 47487770D 10 
0 . 1114 504 15 081777 20D 10 
-0.39179754649029220D 09 
0. 607 817 5 896694 3 HOD 08 

1-8 

0.3397 3558713789260D 01 
-0 .705 826 87 36174 165 OD 03 
0. 355 1645715 16 59940D 05 
-0.76127530737318370D 06 
0.86993044124879000D 07 
-0. 5 95815 485 7765 16 50D 08 
0 . 26 18404 392087 3030D 09 
-0. 76 86705 09 367 976 20D 09 
0.15 3 87 9392418 34 26 OD 10 
-0.21076999149043930D 10 
0.19436719741531690D 10 
-0.115 39 107 69915 8290D 10 
0. 39836531833996820D 09 
-0.60781758986943040D 08 



(N - 14) 



I - 9 

-0.43454135994380020D 01 
0.89977 2 88895 6 5 1130D 03 
-0.4480897 5 5 91674030D 05 
0.94395747735251700D 06 
-0. 105 375 3 640804 1470D 08 
0.7 02056 44801013 5 60D 08 
-0.29957197201463970D 09 
0. 854191659976205 20D 09 
-0.16637140741651880D 10 

0 . 22222615089891290 D 10 
-0.20034483623153590D 10 
0. 1165 5 9225 2464085 OD 10 
-0.39522642489514270D 09 
0 .59347 261881679870D 08 

I - 10 

0 . 5 804710012 117 142 OD 01 
-0.11950398764140700D 04 
0.58456906195874110D 05 
-0.11955987721625510D 07 
0 . 12847 342006460540D 08 
-0.82123426856577030D 08 

0 . 3 36 81964462 39 Shod 09 

-0. 926 84 3 405 44 2477 90D 09 
0. 17 50089836 66 5 44 3 OD 10 
-0. 2 27 6022105 5625 940D 10 
0. 2005 4 346 475 4577 40D 10 
-0.11440295227590730D 10 
0.38140616763727160D 09 
-0 . 5 644086 681409 3220D 08 
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ELEMENTS 0? THE EXPLICIT I.'JVSILE 
IV A PH IX WITH DIVISION BY WEICHTS 



(N = 14) 



I = 11 

-0.828534743267186803 01 
0.16o6931o54741400GD 04 

-0.7968678015b054520B 05 
0.154113155599972203 07 
-0.15521o0o435712500D 08 
0.93602951931Hb650D 08 
-0. 3b 57 2 146 3477 4275 OD 09 
0.9 67 39086 8954 3 56 90D 09 
-0.176 896 35 876 d97900D 10 
0.22410783352539520D 10 
-0.19325345179767670L 10 
0.1C828933169120550D 10 
-0.35565600666^448803 09 
0.51968787273639930D 08 

I = 12 

0.131921826178244403 02 
' -0. 26 176 69 8889117 130D 04 
0.11 37612 98417 3 37 90D 06 
-0.1-4863732642391203 07 
0. 178546024869962O03 08 
-0.100452131915341103 09 
0.372871717932001103 09 
-0.948965638832966603 09 
0.168456591390204803 10 
-0.208516141277610903 10 
0.176512424340394o0B 10 
-0.974397646007497703 09 
0.316122567341825303 09 
-0.457247363567773603 08 



( N = 1 4) 

I = 13 

-u.25929439201462b203 02 
0.4/2054317549156003 04 
-0.131^26455000008603 06 
0.219179309431888303 07 
-0.181703146097125203 08 
0. 9581996543 o587990D 08 
-0.340250269073900603 09 
0.838907871971312803 09 
-0.143464165174180603 10 
0.176868116647374403 10 
-0.147653891687943503 10 
0.806141327716665003 09 
-0.259215490252336403 09 
0.372210506661327203 08 

I = 14 

0 . 908123606831007903 02 
-0.582896639084862403 04 
0.141731478399799803 Oo 
-0.181157314896707403 07 
0.140125621972124603 08 
-0.708148602749950203 08 
0.24448O590251506803 09 
-0.591124991647614203 09 
0.101057153602789803 10 
-0.121578708277531203 10 
0.100674323247962803 10 
-0.546149836395131603 09 
0.174720942033971503 09 
-0.249846126996691203 08 
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ELEMENTS OF THE EXPLICIT INVERSE 
MATRIX WITH DIVISION BY WEIGHTS 



N » 1 5 
I - 1; 

0. 626741147 2 40 37670E 00 
-0.149767 326 413 34850D 03 
0.87991702158178660D 04 
-0.22384410147813O90D 06 
0.30907278460131640D 07 
-0.260707 3505 3359^802 08 
0.14398922790202560D 09 
-0.54295844801453790D j 2 
0.14312370236968300D 10 
-0.266154781029410000 10 
0. 3474 5 206 97496 32 10D 10 

- 0 . 3114234958523696 OD 10 
0.182426 13812 67898CD 10 
-0.62870665172053600D 09 
0.96634843254008150D 08 



1 = 2 

-0.95932021867915070D 00 
0. 2 29246 6110739 235 OD 03 
-0.13462431736329^600 05 
0.34227829881890100D 06 
-0.47221560871994830D 07 
0 .39789495277815740D 08 
-0 .21946 5 3716 4 3265 70D 09 
0 .826 2 35 8 3095 2496 40D 09 
- 0 . 2173832422737238 OD 10 

0.40336319906660950D 10 
-0.52525395858817300D 10 
0.46 945 865 5 20265 5 70D 10 
-0.27413061611607320D 10 
0.94143312509001210D 09 
-0. 144140370287215 30D 09 



( N = 1 5) 

1 V 3 

0.12 3 8087145 95 34810D 01 
-0.29580129552387820D 03 
0.17359797991491830D 05 
-O. 44 O 89 I 6703342427 OD 06 

0.60732983cs14250210D 07 
-0 . 5 10714256 382717 80D 08 
0. 2 80987966439 39410 D 09 
-0 . 10546746851127840 D 10 
0.27650329106754510D 10 
-0 . 51099 3815 5 C103200D 10 
0.66235585027415350D 10 
-0.5 88957 961774 5 3 840D 10 
0.3419557 3 387 8896 30D 10 
-0.1167O455028865680D 10 
0 . 17747 3210826 816 20D 09 

1 = 4 

-0.15132903440408250D 01 
0.36143458277718120D 03 
-0. 2 119059194026 8440D 05 
0.53727739053594980D 06 
-0.73832858170019010D 07 
0.61893331012305830D 08 
-0.33921194402664700D 09 
0.1267 3395 9480705 50D 10 
-0. 3 304614218215 4 370D 10 
0.60700974692866130D 10 
-0.78144193438138940D 10 
0.68962211361328180D 10 
-0.39712890466219770D 10 
0.13434429202464070D 10 
-0.20239306917933810D 09 
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SL.* ENTS OF THE EXPLICIT 18 'SHOE 
i/jaTHIX WITH DIVISION B v Util-}' SS 



( N = 15) 



1 = 5 

Q. 181245 6486 3o 149 300 01 
-0 . 432 ub2 106 5176 1690D 03 
0.24331031296391310D 05 
-0. 64067 995 6 98 3369 90D 06 
0.877 35 6790299 3 29 oOD 07 
-O.73215645062549870D 08 
0 . 39 9036 61011349 5 10D 09 
-0.148106S6760755490D 10 
0.38330647941031090D 10 
-0.698095 564132o27 40D 10 
0 . 0903794 ot>162 35 620D 10 
-0. 777 91485 387 3 2 86 10D 10 
0.44322964323493450D 10 
-0.14o27963179240750D 10 
0.220o 344992b -7 9640D 09 

1 = 6 

-j. 216071642 3482o870D 01 
0.51547223690100050D 03 
-0.30115551313185460D 05 
0.75902806615901360D 06 
-U.1O342913716251820D 08 
0.85762 16147 3356650D 08 
-0.46 380045 5 12 426 4 30D 09 
^.170592788286904700 10 
-0.43 7016 5 136 4934^ 90D 10 
0.7 87 05 6 8407 4 366 690D 10 
-0.991b940^53391o700D 10 
0. 85 5 8005 5 6 66 43 147 OD 10 
-0.48135126963171670D 10 
0.15394035471711490D 10 
-0 .233635 3'413 47084 20D 09 



( N = 15) 

1 = 7 

0.256O9505065780910D 01 
-0.61703749015^204400 03 
0.359428395360953100 05 
-0.9013^0357017199900 06 
0 . 1^19o65 675 5 47 37000 08 
-0.10023ol2S92975b700 09 
U.53625232965057010D 09 
-b.l9463799634265430D 10 
0 .49^45 0c6201087030D 10 
-0. 874 2 ob557 47 374130D 10 
0.10656370799125O70D 11 
-0.922764710195435500 10 
0.5 11 35 85 2 446 70 1240D 10 
-0.l6o40995423502350D 10 
0.241194716946680400 09 

1 = 8 



-0.31420O96437500060O 01 
0.74761738281250160D 03 
-0.43373408203125070D 05 
0.10796483007812500D J1 
-0.1446 428 3 61835 9400D 08 
0 . 117 3 6 246 005 8593 80 D 09 
-0.6lo639408l0547000D 09 
0.22110096937500030D 10 
-0.54916100786377040D 10 
0.957839651323243900 10 
-0.11666231979296060D 11 
0 . 9764670345 117 2080D 10 
-0.5 3-35 100375468370D 10 
0.170587614462b9080D 10 
-0. 24 3696 5 9208984 420D 09 
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ELEMENTS of the explicit inverse 
MATRIX WITH DIVISION BY WEIGHTS 



( N - 1 5) 

1*9 

0 . 3^9 310096 22 29B320D 01 
0.9245 969295 977 5470D 03 
0.53278534088889670D 05 
0.13120350817553940D 07 
0 . 17312 394046 906 410D 08 
•0.13791130423589840D 09 
0. 712038777 357 89710D 09 
•0 . 24897282238044170D 10 
0 . 60497726 835 06 8710D 10 
0.10329636183294670D 11 
0.12352321377318810D 11 
•0.10130480499360400D 11 
0 .54290104 3626 49800D 10 
•0 .17126 26 494903 2 910D 10 
0 .24 11947 16 9466 8040D 09 

I « 10 

0. 49721421265 3706 50D 01 
0.11769003011083430D 04 
0.6 711705 985 7607 340D 05 
0.16244931453902490D 07 
•0.20943608396814680D 08 
0.16235745200514600D 09 
0 .81442 3611974 69660D 09 
0.27681423968320320D 10 
0.6 5 49890909 1008840D 10 
0.10915782300378510D 11 
0.12772292374733760D 11 
0.10273936760151980D 11 
0.5 412082 845 676 86 40D 10 
0.16814912316880170D 10 
0.233635 34134708260D 09 



( N - 15 ) 

I -4o.il 

0 . 66366757805783670 D 01 

-0.15618639807458830D 04 
0 . 87490785 27 1139800D 05 
-0 . 20562089908785 320D 07 
0.25527491884002960D 08 
-0 . 1899921S192076590D 09 
0 . 916776 5 485 2783 980D 09 
-O.30O96847588O33260D 10 
0.69095 89944106407 OD 10 
-u . 11220336871660480D 11 
0 . 128403121265 92 36 OD 11 
-0.10134053583039180D 11 
0.52518637322508010D 10 
-0 .160888667 17 5 5065 OD 10 
0.2 208344 99262795 30D 09 

, 1-12 

-0. 94 6920477150285 70D 01 
0.22038877080534270D 04 
-0 . 1192^5 535 50614790D 06 
0.26504614225451120D 07 
-0.30863745176710530D 08 
0.21688274440203220D 09 
-0. 9 9775 745 93109 866 OD 09 
0. 3150898795 69 21740D 10 
-0.70097464706489320D 10 
0.11094760259234770D 11 
-0 . 124 318509814477 20D 11 
0 . 9641776825 389 9 200D 10 
-0.492430037873&3290D 10 
0. 14 900600482642 920D 10 
-0.202 39 306917 93 3 360D 09 
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ELEMENTS OF THE EXPLICIT INVERSE 
JvO MIX WITH DIVISION BY WEIGHTS 



( N = 1 5) 

I = 13 

0.15^747116376727103) 02 
-0.34193093774b423aOD 04 
0.17021473btt9403o90D 06 
-0.3354^2324160300700 07 
0.3556l79859362oo60D o8 
-0.23330092608529340D 09 
0.10202340^6660^3100 10 
-0.31ul35041do674590D 10 
0.67J05-7S094155080D 10 
-0.103652159413701200 11 
0.114046 37 27 96 97 67 OD 11 

- 0.87158079635 305240 D 10 

0.43 98027 986 503o 2 30D 10 
-0.13175794466cc7660D 10 
0. 17747 3 210526 60670D 09 

I = 16 

-0.296 250436 5 37 50330D 02 
0. 616605819665 35 06 OD 04 
-0. 2.. 6467421o5995170D 06 
0.377686539365636701 07 
- 0.362667971357042600 08 
0.223092702400574403) o9 
-0 . 9 33543404^015 2000D 09 

0.274954783391293900 10 
- 0.580429454293298800 10 

0.882119583756752300 10 
-0.957293103193412600 10 
0 .72363984094249240D 10 
-0.36194492311267530D 10 
0.107 6 5 3205 89 30856 OD 10 
-0 .144 14037026719 9 50D 09 



( N = 1 5) 

I = 14 

0.103765013261067700 03 
-0 .76 202104731 309220D 04 
0.21252066208342310D 06 
-0.312763106964454200 07 
0 . 2s01 335 8416 827 970D 08 
-0.16516^779512135700 09 
0.o72032557768883000 09 
-0.19413456172002270D 10 
0.404032721703793100 10 
-0.607594053914125800 10 

0.654056302260258600 10 

-0 . 491286 5726 945 86 50D 10 
0.24445456450144930D 10 
-0.724181153835^29200 09 
0. 96 6 34843 25 3 9615 40D 08 
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